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Abstract 

The Standard Model Higgs potential becomes unstable at large 
field values. After clarifying the issue of gauge dependence of 
the effective potential, we study the cosmological evolution of 
the Higgs field in presence of this instability throughout inflation, 
reheating and the present epoch. We conclude that anti-de Sitter 
patches in which the Higgs field lies at its true vacuum are lethal 
for our universe. From this result, we derive upper bounds on the 
Hubble constant during inflation, which depend on the reheating 
temperature and on the Higgs coupling to the scalar curvature or 
to the infiaton. Finally we study how a speculative fink between 
Higgs meta-stability and consistence of quantum gravity leads 
to a sharp prediction for the Higgs and top masses, which is 
consistent with measured values. 
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1 Introduction 


The measurements of the Higgs-boson and top-quark masses imply the surprising fact 
that, in the context of the Standard Model (SM) with no additional physics, our universe 
lies at the edge between stability and instability of the electroweak vacuum [1] (see [2] 
for earlier analyses). For the present best fit values of the SM parameters, the Higgs 
potential develops an instability well below the Planck scale, but the proximity to the 
stability region insures that the electroweak vacuum lifetime can be exceedingly longer 
than the age of the universe. 

This intriguing result offers a testing ground for phenomena occurring in the early 
universe. Indeed, the presence of a minimum of the SM potential deeper than the elec¬ 
troweak vacuum raises many cosmological issues: how did the Higgs held end up today in 
the false vacuum? Why didn’t the primordial dynamics destabilise the Higgs held? How 
did patches of the universe with large Higgs values evolve in time without swallowing all 
space? Addressing these questions leads to interesting constraints on early-time phenom¬ 
ena and inhationary dynamics. These constraints are the subject of this paper. Several 
aspects about electroweak-vacuum decay from thermal or inhationary Higgs huctuations 
have already been studied in the literature [3~6], but here we give a comprehensive de¬ 
scription of the phenomenon and reach new conclusions. 

The ehects of the thermal bath during the radiation-dominated phase of the universe 
are twofold. On one side, thermal huctuations can trigger nucleation of bubbles that probe 
Higgs-held values beyond the instability barrier. On the other side, thermal corrections to 
the Higgs potential tend to stabilise low held values, creating an ehective barrier. For the 
observed values of the SM parameters, the latter ehect is dominant and thermal corrections 
do not destabilise the electroweak vacuum, even when the reheating temperature is close 
to the Planck scale [3]. 

More subtle is the issue of the Higgs-held huctuations generated during inhation. In the 
case in which the Higgs has no direct coupling to the inhaton and is minimally coupled to 
gravity (and hence is ehectively massless during inhation), the held develops huctuations 
with amplitude proportional to H, the Hubble rate during inhation. These huctuations 
pose a threat to vacuum stability. For values of H smaller than the height of the potential 
barrier, the Higgs held can tunnel into anti-de Sitter (AdS), according to the Coleman-de 
Luccia bubble nucleation process [7]. When H becomes comparable to the barrier height, 
the transition is well described by the Hawking-Moss instanton [8], which corresponds to 
a thermal overcoming of the barrier due to the ehective Gibbons-Hawking temperature 
T = H/2'k [9] associated with the causal horizon of de Sitter (dS) space. However, a 
more convenient way to compute the evolution of the Higgs huctuations during inhation 
is through a stochastic approach based on a Fokker-Planck equation that describes the 
probability to hnd the Higgs held at a given value h and time t [10]. This approach was 
followed in [3, 6] to derive the probability distribution of Higgs patches in the universe. 
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In this paper, we describe the long-wavelength modes of the Higgs field using a Langevin 
equation sourced by a Gaussian random noise that mimics quantum fluctuations during 
inflation. This method has the advantage of bypassing the problem of choosing boundary 
conditions and it is shown to agree with the results from the Fokker-Planck approach with 
appropriate boundary conditions. 

Quantifying the probability for the existence of a patch of the Higgs held in the SM 
vacuum sufficiently large to encompass our observable universe is a subtle issue, which re¬ 
quires an understanding of how AdS bubbles (with large Higgs-held conhgurations) evolve 
in a dS background, during inhation, and in a Minkowski background, after inhation. As 
correctly pointed out in [6], patches in which the Higgs probes held values beyond the bar¬ 
rier do not necessarily end up in the AdS vacuum, as long as their evolution is driven by 
the stochastic quantum term. Only when classical evolution takes over, the held falls into 
its deep minimum. In [6] it was assumed that these AdS patches rapidly evolve into relic 
defects that are not necessarily dangerous, hence arguing that large Higgs huctuations do 
not pose a cosmological threat. In our analysis, we reach opposite conclusions. 

The evolution of AdS bubbles in an inhationary dS background depends on their size, 
internal energy, surface tension, and initial wall velocity. Depending on the characteristics 
of the bubbles, we hnd a variety of possible evolutions. Bubbles shrink, if they start with 
small radius and low velocity; expand but remain hidden inside the Schwarzschild horizon, 
if the gravitational self-energy of their surface overwhelms the diherence between the 
vacuum energy in the exterior and interior; and expand at the expense of exterior space, 
otherwise. The seemingly paradoxical situation of an expanding bubble of crunching AdS 
space is resolved by understanding the diherence in space-time coordinates on the two 
sides of the wall. While an observer inside the bubbles will experience space contracting 
because of the negative cosmological constant, an external observer will see the surface 
of large bubbles expand. Although we expect that the process of inhation with large H 
will generate a distribution of AdS expanding bubbles, we conclude that such bubbles will 
never take over all dS space. The inhationary space expansion always beats the causal 
expansion of bubbles, efficiently diluting them. 

At this stage, it may seem that the remnant AdS bubbles can be compatible with 
the presently observed universe. The problem starts when we consider post-inhationary 
evolution of the AdS patches in a hat background. The bubble wall keeps on expanding 
at the speed of light and an AdS patch eventually engulfs all space. This means that 
a necessary requirement for our present universe to exist is that the probability to hnd 
an expanding AdS bubble in our past light-cone must be negligible. Unfortunately, we 
cannot make hrm statements about the formation of expanding AdS Higgs bubbles during 
inhation because the answer depends on energy considerations based on the Higgs poten¬ 
tial in the Planckian region. However, our study suggests that it would be very difficult to 
imagine a situation in which all large-held Higgs patches shrink and none expands. There¬ 
fore, barring the presence of AdS Higgs bubbles in our past light-cone is a well-justihed 
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requirement for a viable cosmology. This line of reasoning leads to an interesting bound 
on H, the Hubble constant during inflation, which we compute not only in the case of a 
minimally-coupled Higgs, but also in the presence of a gravitational interaction between 
the Higgs bilinear and the scalar curvature. 

Having established the dangers of patches in which the Higgs held falls into the trans- 
Planckian region, we consider the fate of patches in which, at the end of inhation, the 
Higgs held has huctuated beyond the potential barrier, but has not yet experienced the 
classical evolution that wants to drive it towards very large values. The eventual fate of 
such bubbles is determined by the subsequent thermal evolution of the universe. Thermal 
ehects can rescue such patches of the universe by ehectively pushing the potential barrier 
to larger held values, allowing the Higgs held to relax into its SM vacuum. We study this 
phenomenon during the preheating and reheating stages of the universe, when the energy 
stored in the inhaton oscillating around its minimum is released into thermal energy 
carried by SM particles. In this way we can express the constraint on if as a function of 
the reheating temperature after inhation. 

Finally, in a more speculative vein, we explore the consequences of a conjecture put 
forward in the context of quantum theories of gravity. It has been argued that no for¬ 
mulation of quantum mechanics in dS spaces can be consistent. On the other hand, we 
observe today a positive cosmological constant. The resolution of this conhict between a 
conceptual obstruction and an empirical fact can be found by assuming that the asymp¬ 
totic state of our universe is not dS. In other words, we are only living in a transitory 
situation and today’s dS space will soon terminate. Of course, there are many ways in 
which the universe could escape the allegedly dreadful dS condition, but it is tempting to 
speculate that the instability of the electroweak vacuum is the emergency exit chosen by 
nature. If we take this hypothesis seriously, we obtain a rather precise prediction for a 
combination of the Higgs and top masses in the SM, in good agreement with experimental 
measurements. The predicted strip in parameter space can be narrowed further by the 
hypothesis that the universe must have been sufficiently hot in the past (for instance, for 
allowing some high-temperature mechanism of baryogenesis). 

The paper is organised as follows. In section 2 we address the preliminary technical 
issue of the gauge dependence of the effective potential. The generation and evolution of 
Higgs fluctuations during inflation is studied in section 3, while the subsequent evolution 
after inflation is the subject of section 4. Our speculations on the quantum-gravity pre¬ 
diction of the Higgs mass are discussed in section 5, and our results are summarised in 
section 6. The details of the general-relativity calculation of the evolution of AdS bubbles 
in dS or Minkowski backgrounds are contained in the appendix. 
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2 Gauge dependence of the SM effective potential 


The critical Higgs mass below which the SM Higgs potential becomes unstable is gauge- 
independent; however the instability scale of the SM potential, e.g. hmax, defined as the 
Higgs field value at which Ves{h) is maximal, is gauge-dependent (as recently emphasised 
in [11,12]). Therefore, one has to be cautious in extracting from the potential a physically 
meaningful scale associated to the instability. 

There is a number of ways in which one can try to identify scales that track the poten¬ 
tial instability and are gauge-invariant because expressed in terms of extrema.^ However, 
our cosmological computations will employ the full SM effective potential also away from 
its extrema, so we are confronted with the issue of the gauge-dependence of the effective 
potential shape, an old topic much debated in the literature.^ In dealing with this issue we 
follow a pragmatic approach. First, we insist on calculating physical quantities, that can 
be proven to be gauge independent. Second, we make sure that the approximations we 
use in those calculations are consistent, in the sense that any residual gauge dependence 
is smaller than the precision of our approximations. 

The gauge-independence of our results is ultimately based on the Nielsen identity that 
describes how the effective action depends on the gauge-fixing parameters and how to 
extract out of it gauge independent quantities. Let us briefly discuss how this works. 

The fact that the Higgs effective potential V{h) depends on the gauge parameters 
(generically denoted as follows from the fact that the effective action S^s itself is a 
gauge-dependent object. In spite of this, as is well known, both the potential and the 
effective action are extremely useful and physical quantities extracted from them (like 
particle masses, S'-matrix elements, the vacuum energy density, tunnelling rates in the 
case of metastable vacua, etc.) turn out to be gauge-independent, as they should. 

^For example, this could be done through the scale of a higher-dimensional operator (with n > 4) 
that, added to the SM, cures its instability. 

^We summarise some of the main approaches here. Nielsen [13] proved that the gauge dependence of 
the effective potential can be reabsorbed by a re-definition of the fields. Tye [14] found that the effective 
potential and the effective kinetic terms are separately gauge-invariant if the perturbative expansion is 
performed by decomposing the Higgs doublet into the physical Higgs field h and the 3 angular coordinates 
TT, such that the Goldstone fields tt are massless at any value of h, not only at the extrema of the potential. 
The results then agree with the unitary gauge. Buchmuller et al. [15] computed the effective potential 
in terms of the gauge-invariant combination claiming a gauge-invariant effective potential; again 

this selects the radial mode of the Higgs doublet such that Goldstone are always massless. Schwartz et 
al. [12] argue that finding a gauge-invariant definition of the effective potential is a misguided enterprise 
and that the contribution of Goldstone bosons should be neglected (at leading order) in a consistent 
perturbative expansion around the field value at which A = 0. 
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2.1 Gauge (in)dependence of the effective action 


The Nielsen identities [16] tell us that the gauge dependence of the effective action can 
be compensated by a local held redehnition. In other words, different gauges describe the 
same physics in terms of different coordinates in held space (leading to diherent potentials 
but also to diherent kinetic terms). Particularising to cases with Higgs background only, 
one has 

( 1 ) 

where K[h{xy\ is a functional of h that can be found in [16]. 

One immediate consequence of the Nielsen identity is that the action evaluated on a 
solution of the equation of motion for h, SSes/Sh = 0, is gauge-independent. We also 
see that the gauge-independence of the extremal values of the ehective potential follows 
directly by applying the previous general fact to constant held conhgurations. 

Writing the ehective action in a derivative expansion 



-v{h) + ^z{h){d,hY + o{d^) 


( 2 ) 


we can hnd a series of Nielsen identities for the coefficients of this expansion [17,18] by 
expanding in the same way K[h] and 6Ses/Sh in (1), as 


K[h] = C{h) + D{h){d^,hf - d^^[D{h)d^,h] + 0 ( 9 ^ 


5S, 


eff 


5h 


= -V + -Z'{d,hf - d^[z{h)d,h] + o{d ^), 


(3) 

(4) 


where primes denote h-derivatives. The higher-order derivative terms are expected to 
be suppressed by an energy scale, which can be as low as the value of the Higgs held, 
and by a one-loop factor, which is the same for any order in derivatives. Therefore, our 
derivative expansion is valid only when the gradient of the Higgs held is smaller than the 
homogeneous value of the held under consideration. 

At the lowest order in the derivative expansion, we hnd the expression for the gauge- 
dependence of the ehective potential 


dV 

e^ + c'(h)w = o. 


(5) 


which, as anticipated, ensures the gauge-independence of the values of the potential at its 
extremal points. This Nielsen identity also tells us that the explicit .^-dependence of the 
potential can be compensated by an implicit ^-dependence of the held as: 

dh 




so that dV/d^ = 0. At order 0{d‘^) we get 

dZ 


= -CZ' - 2ZC' + 2DV' + 2DV" 


( 6 ) 


(7) 
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where we suppressed the h dependence of all functions. 

It is useful to consider the order in weak gauge couplings (denoted generically by g in 
this paragraph) of the different functions that appear in the previous identities [17]. As we 
will be interested in the potential region where the Higgs quartic coupling gets negative, 
eventually inducing a new minimum radiatively, we will use the counting A ~ g'^. The 
function C{h) starts at one loop and is 0{g‘^). The Nielsen identity (5) then implies 
that the ^ dependence of V starts at O(g^). On the other hand, the Nielsen identity (7) 
implies that the ^ dependence of Z starts at O(g^), with the terms involving D and D 
being of higher order in g. As we will see in the next subsection, it will be sufficient for 
our purposes to deal with the dominant ^ dependence of the potential so that we will 
neglect the effect of the subleading D and D terms in what follows, as in [17]. 

Let us next consider the ^-dependence of the equation of motion for h which we write 
using 

“ 

and its solutions h{x), which satisfy EoM[h] = 0. It is straightforward to show that 

_ = 0 , (9) 

h=h 

up to 0(5"^) corrections, provided we use ^dh/d^ = C. In principle, one can continue 
the check of the gauge invariance of the equations of motion iteratively up to inhnite 
order in the number of derivatives.^ This means that, if some h^{x) solves the equation of 
motion for some choice of ^ and we shift ^ dfy the shifted solution is h^{x) -|- dh^{x) 

with dh^{x) = C(h^)d^/^. In other words, the held rescaling that can balance the effect of 
changing ^ in the effective potential is the same held rescaling that applies to the solutions 
of the equation of motion for diherent fy 


? ^EoMlftl 



The same rescaling works for the Fokker-Planck and Langevin equations that we will 
use later on to describe the Higgs huctuations during inhation. These equations take the 
form, LANGEVlN[hL] = 0 and FoKKERPLANCK[P(h, f)] = 0, with 


LANGEV.Nlft|^^/^f 


( 10 ) 


and 


d (H^ P \ 1 PH'l] 1 dP 

^ x^Vz) ^ 3pWJ J ~ Vz^ ■ 

( 11 ) 

^Using the previous identities one can also check, to all orders, the ^independence of the scalar 
physical mass = V "evaluated at the minimum of the potential, as indicated. 


1 d 1 


FOKKERPLANOK|P(;,,i)l^^^|^ 














Here P{h, t) dh is the probability for finding the Higgs held in the inhnitesimal interval 
between h and h + dh at time t during inhation. The fact that P is a probability density 
explains why P enters in eq. (11) through the ratio P/a/Z. 

Using these expressions and the ^-dependence of U(h), Z{h), and h as described in 
eqs. (5)-(7), we get 


(where hi and P are solutions of the Langevin and Fokker-Planck equations, respectively) 
up to corrections that can be shown to be subleading.This shows once again that if we 
have a solution of the Langevin equation for a given value of we automatically obtain 
a solution for ^ -l- by the shift h{^) -|- C{h{^))d^/^. For the Fokker-Planck equation, a 
solution for general Z is formally related to a solution for Z = 1 again by a held rescaling, 
with P(h, t)/^Jz{h) = Pc{hc{h)), where the relation between h and the canonical held he 

follows from dhc/dh = ^Z{h). As a result, the integrated probability is independent of 
the held rescaling: 


^ Langevin [h] 


h=h T 


= 0 , ^FOKKERPLANCK[P(h, f)] 


= 0 , 


( 12 ) 


P=P 


Pe{he)dhe = Pe{he{h)) ^^)dh = P{h)dh . (13) 

This implies that the probability of hnding the held beyond hmax after a given number 
of e-folds is a gauge invariant quantity: although the value of hmax depends on the ^- 
change of the ratio P{h,t)/Z{h) corresponds to the same held-rescaling h{^) — h{^) + 
C{h{^))d^/^ and leaves the integrated probabilities unchanged. 

To sum up, the key idea is that a change in a given gauge parameter ^ is equivalent to a 
redehnition of the Higgs held, which should leave physics invariant. The ehective potential, 
the equations of motion for the Higgs held and the Fokker-Planck and Langevin equations 
enjoy a sort of “covariance” under changes of the gauge parameters. The equations are 

'^Here we are explicitly using the derivative expansion previously introduced to derive the gauge 
transformation properties of F, Z and h. Indeed, the Langevin and Fokker-Planck formalism represent 
a truncation of the theory at the lowest order in derivatives, where the approximation is justified by the 
smallness of the gradient of the field with respect to the Hubble parameter. In the rest of the paper, we 
will be using these equations to describe evolutions of the Higgs field for values of the Hubble parameters 
even quite larger than the Higgs vev itself. For this reason, we cannot naively apply the zeroth order 
truncation in derivatives of the effective action, because, as we discussed, the derivative expansion is 
suppressed only by the Higgs vev. However, since the higher derivative corrections are suppressed by 
at least a one-loop factor and are not log-enhanced (at one-loop), a consistent truncation is to use the 
Langevin and Fokker-Planck equations as derived from an effective action where the only corrections 
that are included are the non-derivative, leading log-enhanced, ones. This will be how we will use the 
Langevin and Fokker-Planck equations in the rest of the paper. 

^ Note also that, concerning the dependence on the renormalisation scale p, one can show the p- 
independence of EoM[/i], Langevin[/i] and FokkerPlanck[P(/i, t)] just making use of dV/dp = 0, 
dh/dlogp = 'yh, d{dV/dh)/dlog p = —j{dV/dh), dZ/d\ogp = —2'yZ, etc. 
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changed in such a way that the change induced in their solutions is just a common held 
redehnition dictated by the Nielsen identity. 

2.2 Effective potential including only log-enhanced corrections 

For the previous appealing properties to hold, the interplay between the effective potential 
and the kinetic term in the effective action is crucial. For the SM case at very large held 
values we write 

= + (14) 

where the ellipsis denotes higher derivative terms and both Z{h,^) and Xes{h,^) include 
radiative corrections and depend on As usual, it proves convenient to use a canoni¬ 
cally normalised Higgs field hcan(h,0 dhca.n/dh = and to re-express the effective 
Lagrangian in terms of h-cam obtaining 

(d h 

=^efr = - A,an(hcan, 0^ + ' ' ' (15) 

In terms of the canonical field all the equations become simpler as we do not have to 
drag the Z factor around. An additional bonus is that the residual ^ dependence in our 
approximations will be significantly reduced. Let us see how this works examining the 
gauge dependence of the effective potential. The coloured dashed curves in fig. 1 show 
Aefr(h, ^), which tracks the large field behaviour of the SM effective potential, as computed 
at next-to-leading order (NLO) accuracy in the Fermi ^ gauges [11]® further improved by 
performing a resummation of IR-divergent Goldstone loops [19]^. We take into account 
the running of and ^2 (gauge-fixing parameters for hypercharge and SU(2)i), assuming 
a common value ^ renormalised at Mt. We confirm that Aefr(h,.^) significantly depends on 
the gauge-parameter The black curves in fig. 1 show Acan(hcan,0) again computed at 
NLO in Fermi ^ gauges: we see that the dependence on ^ almost completely disappeared 
— all black curves almost merged into a single curve. 

One can explain analytically why the gauge dependence approximately cancels out by 
looking at the dominant corrections enhanced by large logarithms, which are resumed by 
solving the RG equations and setting the RG scale ft around the field value of interest: 

Aeff(h, e) ^ ^ h), Z,s{h, 0 ~ ^ (ig) 

where r = /i^.7 d In /i, 7 is the gauge-dependent anomalous dimension of the Higgs field, 
and A(/ 2 ) = /3xdlnj2 is the running quartic coupling, where /3 a is gauge-independent. 

® At this level of approximation the potential has a residual dependence on the RG scale p, comparable 
to the gauge dependence. We here adopted the choice p = he^ that minimises the error. 

^Around the minimum of the potential, this resummation becomes equivalent to the expansion of [12], 
which makes the truncation of the potential compatible with the Nielsen identity. 
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Figure 1: The dashed curves show the effective quartic coupling (left) and effective SM 
potential (right) computed at next-to-leading order in a generic Fermi -gauge. The thick 
red dashed curve corresponds to the Landau gauge, f = 0. The right handed panel shows 
that the height of the potential barrier is only approximately gauge-independent (a measure 
of the residual gauge dependence). The black continuous curves show the same potential 
expressed in terms of the canonical field h r»n ■' the gauge dependence in the potential gets 
compensated by the gauge-dependence of the kinetic term, such that the continuous curves 
nearly overlap. 

In this leading-order (LO) approximation one has 

hcan ~ Acan(hcan) ~ Kfi ~ > (17) 

which is gauge-independent because the RGE for A and all other couplings of the theory 
are gauge-independent. The order-of-magnitude gauge dependence of hmax found in [11] 
disappears because it is almost entirely due to the RG factor F. 

Such LO cancellation has been noticed before, see e.g. [20,21], Note however that 
the field redefinition dictated by the Nielsen identity we discussed earlier and the field 
redefinition required to make the field canonical are the same only at LO. Moreover, 
the field redefinition from the Nielsen identity becomes considerably more complicated 
at NLO [13]. Its use to define a “gauge-independent” potential is in fact equivalent to 
choosing a particular gauge and therefore does not solve the problem of how to extract 
gauge-invariant quantities out of the effective action. For this reason we refrain from 
attempting to use it as a way of dehning a gauge-invariant potential and simply use the 
canonical held as a way of reducing the residual gauge dependence of our results. 

The previous discussion has been carried out in Fermi gauge at NLO to help us clarify 
the issues related to gauge invariance. Having understood them, we can now use the 
state-of-the-art computation of the effective potential in the Landau gauge = 0) with 
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NNLO accuracy (2 loop finite corrections and 3 loop RGE corrections) [1] combined with 
the use of a held redehnition to make the held canonical (taking into account the ehect 
of potentially large logarithms in Z). Using that canonically normalised Higgs held h, 
in the region around the top of the barrier, the SM Higgs potential can be analytically 
approximated as 


Heff(h) -6 In 



T ’ 


(18) 


where hmax is the held value at which Kfr(h) takes its maximal gauge-invariant value 
Kff(h ma x) = Using the value b ^ 0.16/(47r)^ for the {3 function of A around hmax, 

we hnd hmax = 5 x 10^° GeV for the present best-ht values of Mt, Mh and a^. Although 
this value of hmax is computed in Landau gauge and it would be slightly diherent in other 
gauges, the reader should keep in mind that the results we present in the following sections 
are gauge-invariant even if for convenience we express them in terms of hmax-* 


3 Higgs fluctuations during inflation 


The instability of the Higgs potential leads to an interesting dynamics during inhation. 
We focus on the relevant radial mode h = of the Higgs doublet. If the Hubble 

constant H is large enough, h huctuates beyond the potential barrier. If the true vacuum is 
deep enough, inhation stops in the regions where the Higgs falls, while inhation continues 
in the (possibly rare) regions where accidentally h < hmax- In this section we compute 
the probability of the possible outcomes at the end of inhation, while in the next section 
we will discuss what happens after inhation. 

In the absence of a large Higgs mass term, the evolution of the long wavelength modes 
of the h held is controlled by the Langevin equation [22] 


dh 1 dV{h) 
dt 3H dh 




where rj is a. Gaussian random noise with 




dvr 2 




(19) 


( 20 ) 


It is important to realise that eq. (19) is valid only if the positive ehective mass squared 
V"{h) of the Higgs held is light enough compared to H^. Only under these circumstances 
the long wavelength super-Hubble huctuations of the Higgs held are generated. On the 

®Alternatively, we could choose other scales associated (more indirectly) to the instability which are 
explicitly gauge invariant. One could be the renormalisation scale /io at which the quartic Higgs coupling 
A crosses zero; another choice is the scale at which the one-loop radiatively corrected Goldstone 
mass is zero (as used in [12]). For the same central values above we get /tq = 1 x 10^° GeV and 
^jix = 4.5 X IQio GeV. 
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contrary, if V''{h) > 9iJ^/4, the resulting power spectrum of Higgs fluctuations is both 
suppressed by exp(—21/"(h)/if^) and by the fact that the spectrum is strongly tilted on 
the blue side [23]®. 


3.1 Higgs fluctuations during inflation for = 0 

It is convenient to rewrite eq. (19) replacing time t with the number of e-folds N = Ht, 
and to normalise the Higgs held and its potential in units of the Higgs value h = hmax at 
which V (hjnax) = hmax IS maximal, 


h 


h 


and 


vCh) 


V 

'^max 


—6 In 



T‘ 


After these redehnitions, the Langevin equation in eq. (19) becomes 


dh dVjh) 

dN dh 


fj{N) 


where the noise 'r]{N) obeys 


{r!{N)ri{N')) = 


H 


27lhr. 


2 


5{N-N'). 


( 21 ) 


( 22 ) 


(23) 


One can now numerically generate random realisations of the Higgs evolution in N, 
in steps of dN, as 

h{N+ dN) = h{N)-^V'{h)dN+ r (24) 

where r are random numbers extracted from a Gaussian distribution with zero mean and 
standard deviation a = Hy/dN / (27rhjnax)- 

Indeed, for h ^ a, the same result is reproduced by the analytic solution to the 
Fokker-Planck equation for the probability P{h, N) of hnding the Higgs held at the value 
h after N e-folds of inhation. 


dP 


\ 

d , 

^ V' 

M ~ 


1 

1 00 

dh ' 



taking V = 0 and boundary conditions at h = ±cx). 

Figure 2 shows the resulting probability density of the value of h after = 60 e- 
foldings, starting from h = 0 at the beginning of inhation. The result, a quasi-Gaussian 
distribution, has a simple interpretation. Given that the quartic Higgs coupling vanishes 
around the instability scale, for a large range of Higgs values around the instability scale 
the classical evolution (sourced by the gradient of the potential) is negligible with respect 

^Indeed, the solution of the Klein-Gordon equation for a spin 0 particle with mass m in de Sitter goes 
like exp(—7r/i/2)(—fcr), where fj, = y/m'^/H'^ — 9/4 and r is the conformal time. 
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H = /?max, jO(l^l>^max) = 0.42, p(\h\^oo) - 0.00016 



ll — /l//^max 


Figure 2: Random distribution of the Higgs field h = /i/Zimax after N = 60 e-folds of 
inflation with Hubble constant egual to the Higgs instability scale, H = /imax- The blue 
dashed curve is the V = 0 Gaussian approximation of eg. (26). The red curve is the SM 
Higgs potential V(h), in arbitrary units. 


to the quantum evolution (sourced by the random noise t]). As a consequence, even 
assuming that the Higgs starts from h = 0, the held h acquires a Gaussian distribution 
with zero mean and variance that grows with N\ 


P{h,N) 



(26) 


The distribution shown in £g. 2 maintains its quasi-Gaussian shape also for values of 
the Higgs held well above hmax- Therefore, during inhation, the Higgs held can huctuate 
above the barrier without being sucked into the negative-energy (AdS) true vacuum. The 
only regions where the Higgs falls into the true minimum are those where h huctuates to 
held values so large that the potential slope can no longer be neglected. 


Regions that fluctuate above the potential barrier 

The minimal probability that the Higgs ends up beyond the top of its potential barrier 
after N e-folds is 

p{\h\ > hmax) ~ 1 -erf(^^^^^^y (27) 

This probability, obtained by integrating the Gaussian distribution for \h\ > h-max, is 
minimal because it corresponds to the initial condition h = 0. Shifting the peak of the 
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yv = 60 e-folds, -0 



Figure 3: Minimal probability that, after N = 60 e-folds of inflation, the Higgs fluctu¬ 
ated above the SM potential barrier (orange curve), or fall down to the true minimum 
(red curve). The continuous curves are the numerical results; the dashed curves are the 
analytical approximations presented in the text. 


distribution to a non-vanishing value of h will only increase p{\h\ > hmax)- The solid 
orange curve in £g. 3 shows our numerical result for this probability as a function of the 
Hubble constant during inflation, in units of the Higgs instability scale if/hmax- The 
dashed orange curve corresponds to the analytic expression in eq. (27), which is evidently 
an excellent approximation. 

No constraints arise if, after inflation, the regions with \h\ > hmax fall back to the SM 
minimum, pushed by thermal effects (see section 4). If instead, after inflation, the regions 
with \h\ > h-max fall down into the true AdS minimum, then their probability should be 
smaller than so that it is unlikely to find the Higgs away from its EW vacuum in 

any of the ~ causally independent regions that are formed during inflation and that 
constitute the observable universe today. Using 1 — erf(a;) ~ /y/Ax for large x, this 

condition implies 



(28) 


Regions that fall to the true minimum during inflation 

The approximation of neglecting the scalar potential V, which led to the quasi-Gaussian 
distribution of the Higgs held values, breaks down at large h. There, the gradient of the 
potential dominates over quantum fluctuations, and h falls down to its true minimum 
already during inflation. The solid red curve in fig. 3 shows our numerical result for such 
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probability. 

We can analytically estimate the probability for h to fall into its true vacuum after N 
e-folds of inflation. We first consider a potential V = Ah^/4 with constant A and assume 
that the bulk of the Higgs held probability distribution is still given by the Gaussian in 
eq. (26), cut at large held values. The location of the cut is estimated by demanding that 
the classical evolution becomes more important than the quantum huctuations [6]. This 
can be quantihed by requiring that the second term in the right-hand side of eq. (25) 
dominates over the hrst one, 


d / V 



> k 


dh^ 



(29) 


where A; is a fudge factor and P is given in eq. (26). Equation (29) implies that the 
Gaussian distribution must be cut for > 3fcW^/(2|A|iV), and values of h that satisfy 
this inequality are sucked into the true minimum. 

Therefore, the probability of falling to inhnity is exponentially suppressed for small 


|A|: 


p{\h\ —)■ cx)) 1 — erf 



(30) 


Such probability satishes p{\h\ —>■ cx)) < for |A| < kn'^/N^. 

Gonsidering now the more realistic case of the SM potential with a running coupling 
X{h) = —b\n{h‘^/h‘^^^^/e), we hnd 


p{\h\ —)■ oo) cs 1 — erf 


/ tta/SA: a 

VivTH/ ’ 


where 


B = FL 


/ 3kH^ 

1 ^ 26 ML 


(31) 


where PL is the ProductLog function. This analytic approximation of p(|h| —)■ cx)) is 
shown in Eg. 3 as the dashed red line and agrees well with the numerical computation, 
once we £t the fudge factor to be fc = 2.6. Next, we need to extrapolate the analytic 
approximation to probabilities much smaller than those that can be computed numerically. 
The probability p{\h\ —)■ cx)) is smaller than for 


H 


71 




^k/2bN^ 


^ 0.045 


(32) 


3.2 Higgs fluctuations during inflation for 0 

Higgs fluctuations during inflation can get damped if the Higgs doublet during inflation 
acquires an effective mass m. Various effects can contribute to such mass: 

1 . a quartic term in the potential, which describes a coupling between 

the Higgs and the inffaton 0, generates during inflation an extra contribution 
to the Higgs mass; 
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Figure 4: Running of the Higgs coupling to gravity os a function of the renormalisation 
scale in the SM, for different initial conditions at the Planck scale. The dashed horizontal 
lines correspond to the special values fn = —1/6 and = 0. 


2. a decay of the inflaton into SM particles can generate a non-vanishing temperature 
during inflation. Such decay are kinematically blocked when SM particles acquire 
a thermal mass larger than the inflaton mass, of order H. Thereby the Higgs can 
acquire a mass 

3. a non-minimal Higgs coupling to gravity, —fnl^nl^R contributes as = fnP = 
-U^hHT 


These contributions to would have qualitatively similar effects. In our quantitative 
analysis we focus on the latter effect because the presence of the fn term is unavoidable: 
even if = 0 at some energy scale, SM quantum corrections generate a non-vanishing 
value of fn at any other energy scale. Indeed, ignoring gravity, the one-loop running of 
is given by 


df,H _ 'CiT + 1/6 
din ft m? din ft 


in + 1/6 

(47r)^ 




9 

10 


Qi -I- 12A_h' 


+ ■ ■ ■ 


(33) 


where ft is the renormalisation scale. The RGE for the Higgs mass parameter is known 
up to 3 loops in the MS scheme (as summarised in [1]), and we have shown here only 
the leading term. The SM couplings are such that din |m^|/dln/i is positive (negative) 
at energy roughy below (above) 10^° GeV. The evolution of fu for different boundary 
conditions at Mpi is shown in Eg. 4; it has a fixed point at the conformal value = —1/6. 
Notice that this value is not special for our analysis because it does not recover conformal 
invariance, which is broken at the level of the SM Higgs effective potential. 
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H = /Zmax, P(l^l>/lmax) = 0.075, p(\h\->oo) = 0. 



It — 

Figure 5: Random distribution of the Higgs field h = /i/Zimax after N = 60 e-folds of 
inflation with Hubble constant egual to the Higgs instability scale, H = /imax? and for 
fn = —0.01. The blue dashed line is the Gaussian approximation of eg. (37). The red 
curve is the Higgs potential Vsm(^) — in arbitrary units. 


We consider the following action 


S = 









-V + ■■■ 


(34) 


where V ~ V((/)) + is the scalar potential of the Higgs and of the inflation and 

Mpi is the reduced Planck mass. We use the approximation that, during inflation, the 
inflaton potential is constant H(0) ~ 

The fu coupling of the Higgs to gravity affects the scalar potential during inflation by 
inducing an effective Higgs mass term m? = fnR = which can 

stabilise the Higgs potential and suppress Higgs fluctuations. As explained after eq. (19), 
Higgs fluctuations are damped if fn < —3/16. For —3/16 < fu < Higgs fluctuations 
are still present, but become less dangerous than in the case of vanishing 

course, one could also envisage other operators coupling the Higgs field with gravity, e.g. 
pi?^/Mp[. However, in most models of inflation, the Hubble parameter squared decreases linearly 
with the number of e-folds W till the end of inflation. Therefore, for ^ (Hj/Mp)){Ng/Nj), where 
Hj is the initial value of the Hubble rate when inflation starts and Nj is the total number of e-folds, the 
higher-order operator becomes negligible. This condition becomes easier and easier to satisfy as inflation 
proceeds. 
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Figure 6: ^45 a function of fu o,nd the Hubble constant in units of the instability scale /imax 
(and for N = 60 e-folds of inflation), we show the three regions where: the probability for 
the Higgs field to end up in the negative-energy true minimum is larger than e~^^ (red); 
the probability for the Higgs field to fluctuate beyond the potential barrier is larger than 
g-3Af )orange); the latter probability is smaller than e~^^ (green). Higgs fluctuations are 
damped for < —3/16. The uncertainty on the orange/red boundary corresponds to a 
fudge factor 1/3 < k < 3. 


Neglecting the small Higgs quartic coupling, adding the effective Higgs mass term 
= —12C,hH‘^, and assuming the ansatz of a Gaussian distribution with variance (h^), 


P{h,N) 




(35) 


the evolution of is obtained from the Fokker-Planck equation (25) and becomes, at 

h^ > (h2). 


?M. = _!!!!!/;,+ iff 

dN iW' ' 4x2 




— exp 


/ 2m^N\ 

3m )' 


(36) 


If imf < 0, the variance grows exponentially with N. If m? > 0, the Higgs probability 
distribution approaches, after a few e-folds, the limiting distribution given by^^ 



H 


(37) 


This is to be compared with eq. (26), which holds for = 0. Figure 5 shows a numerical 
example: already for = —0.01 the variance is significantly reduced. 

^^This is larger than what is obtained by naively assuming a Hawking temperature T = iJ/27r. 
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Using eq. (37), we obtain the following bounds on H from the request that the prob¬ 
abilities of the Higgs fluctuating beyond the barrier (|h| > hmax) or falling into the true 
minimum {\h\ —)■ cx)) are less than e~^^: 



(38) 


(39) 


These bounds are the analogues of eqs. (28) and (32), which are valid for = 0. An 
order-one fudge factor k can be similarly introduced such that eq. (39) closely agrees with 
the numerical result; however k depends on and thereby differs from what we previously 
discussed in the limit = 0. Taking into account how stabilises the potential results 
in a more complicated, but numerically similar, analytic expression. Furthermore these 
approximations needs to be extrapolated down to probabilities smaller than those that 
can be compared to the numerical result. Conservatively estimating the uncertainty by 
varying l/3<fc<3, in fig. 6 we summarise the situation by showing the regions of 
and H/hmax where the bounds in eq.s (38) and (39) are satished for N = 60. 

In the presence of a potential coupling between the Higgs and the 

inflaton 0, during inflation one has an extra contribution to the Higgs mass, 

This term has a similar effect as the inflationary mass discussed above. The Higgs h has 
no inflationary fluctuations as long as m > 3H/2. However, in general changes during 
inflation in a model-dependent way. Considering, for example, large-held inhation with 
a quadratic potential, one has 0 = 2M'p\\/Ni — N during inhation, where Nj > 60 is the 
total number of e-folds. Inserting w? = 4A/i<^M|i(A'/ — N) into eq. (36) one hnds that the 
maximal Higgs huctuation is achieved at the end of inhation and is Planck suppressed: 



(40) 


3.3 Bubble evolution in de Sitter spacetime 

After having computed the probability for inhationary huctuations to form regions where 
the Higgs held lies at its true minimum, the next question we have to address is how these 
regions evolve. In the literature one hnds conhicting statements about the evolution 
of AdS regions in an inhationary background. One point of view, based on hat-space 
intuition, is that AdS regions should expand because their interior has lower energy than 
the exterior. A diherent point of view is that, since AdS space eventually contracts, 
regions in which the Higgs lies at its true minimum will shrink, possibly leaving some 
almost point-like relics, which are nevertheless efficiently diluted, and thus made harmless, 
by the inhationary expansion of space. We will show that addressing the question about 
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the fate of AdS regions involves a nnmber of non-trivial and connter-intnitive issnes raised 
by general relativity. 

First, gravitational energy contributes to the total energy budget. Second, an AdS 
region might expand, while remaining hidden behind a black-hole horizon. Third, the 
interior AdS space is dynamically unstable [7]: when described in cosmological FRW 
coordinates it reaches a ‘big-crunch’ singularity in a hnite amount of internal time of 
order where —Idn < 0 is the internal cosmological constant. If space is empty, 

this is just a coordinate singularity (AdS can be continued using better coordinates); if 
space is hlled by a background held (for example the Higgs held), its huctuations grow 
until the energy density becomes inhnite and a physical singularity appears. Furthermore, 
the AdS geometry has a timelike boundary, such that the evolution cannot be predicted 
after the bubble wall reaches the boundary, unless additional boundary conditions are 
imposed there (in other words, information must how in from inhnity). As a result, a 
Cauchy horizon appears in the interior of the bubble. It is expected that, within the full 
theory beyond the thin-wall limit, a physical spacelike singularity must develop before 
the Cauchy horizon [24]. This conhrms the expectation that the AdS bubble is unstable. 

In order to clarify all these issues we performed a careful (and somewhat lengthy) 
general-relativistic computation, described in appendix A. Here we summarise the main 
points. 

In order to make the problem tractable analytically, we assume a spherical AdS region 
(that we thereby call ‘bubble’), separated from the outside space by a thin wall with 
constant surface tension a. The matching of the external and internal geometries requires 
the presence of such a wall with nonzero energy density. The fate of the AdS bubble is 
then determined by computing the motion of the wall separating the AdS interior from the 
external space (de Sitter during inflation and Minkowski after inflation). In the thin-wall 
approximation, the motion of the wall is determined by junction conditions that relate 
the extrinsic curvature on each of its sides [25]. The bubbles that we now compute are 
more general than those that arise from vacuum decay with zero total energy, already 
studied in [7]. The basic elements of our calculation are the following. 

1 . The space inside the bubble is assumed to be an empty spherical region of AdS 
space with metric 

ds^ = -fin{r) drj"^ + —— r < R, (41) 

hn{r) 

expressed in global coordinates. Here fin{r) = 1 -|- is the usual AdS solution, 

with vacuum energy —Idn corresponding to the length scale 1/^fn = SvrGldn/S. 

2. The space outside the bubble is described by the metric 

ds^ = -fontir) dt^ ^ r > R, (42) 

Jout{r) 
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where fout{r) = 1 — ~ 2GM/r describes a Schwarzschild-de Sitter (SdS) 

spacetime, with G = l/(87rMp[). Here M is the mass of the bubble as seen by 
an outside observer, living in an asymptotically de Sitter space described by the 
length scale l/-^out = 87rG14ut/3 = As discussed later, the metric in eq. (42) 
also describes the case of the asymptotically flat spacetime produced after inflation, 
which is obtained in the limit Kut —^ 0, so that /out(’") = 1 — 2GM/r. 

Note that, for £out 3> GM, the SdS spacetime contains two horizons, corresponding 
to the zeros of the inner (Schwarzschild) horizon at r 2GM and the outer 

(de Sitter) horizon atr ^ £out- The corresponding Penrose diagram is depicted in the 
right panel of £g. 7. It is a combination of the diagrams for the Schwarzschild and 
de Sitter spacetimes [9]. Thick blue lines denote curvature singularities, the dashed 
lines horizons and the dotted lines conformal inhnities. The two thin vertical lines at 
the ends of the diagram indicate that the pattern is repeated indefinitely on either 
side. 

3. The two regions are separated by a domain wall with constant surface tension a. 
The metric on the domain wall can be written as 


ds^ = —dr'^ + R‘^{T)dfll, (43) 

where R{t) denotes the location of the wall in both coordinate systems (41) and 


(43) 


(42). The evolution is expressed in terms of the proper time r on the wall. In the 
full problem, a is given by the kinetic and potential energy of the Higgs held, and is 
different for each Higgs conhguration. Within our approximation, all energy stored 
in the Higgs potential goes into the motion of the wall, leaving the AdS interior 
empty. 

The detailed calculation described in appendix A shows that the (naively positive) differ¬ 
ence between the energy in the exterior (iout) and the interior {—i^) of the bubble that 
controls whether the bubble expands or contracts receives a gravitational correction —k^, 
so that the relevant parameter is the quantity: 



(44) 


As discussed after eq. (90), the contribution ~ can be interpreted, from a Newtonian 
point of view, as the gravitational self-energy of the wall. The motion of the wall can be 
described as the Newtonian motion of a point particle 



2 


(45) 


in an effective ‘potential’ given by 



2 


( 46 ) 
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II 


r = 0 r = oo 



Figure 7: Penrose diagram describing an AdS bubble that expands in a dS 
space. The black curve denotes the thin wall separating the two phases; the true space 
is obtained by patching the region in white in the left panel (AdS interior of the bubble) 
with the region in white in the right panel (Schwarzschild-dS exterior of the bubble). The 
dashed lines denote the various horizons, while the thick blue line in the left panel denotes 
the AdS singularity (‘crunch’). 


where R = pR is a rescaled dimensionless coordinate that describes the position of the 
wall as a function of a rescaled dimensionless proper time r = ‘iKrj'y^. The various 
constants are given by 


Anf 3 |A| 

^ ^ 2GM’ 


e = sign A, 



E = 


K 




(47) 


The possible types of bubble evolution are discussed in detail in appendix A. The com¬ 
plete analysis is performed for an asymptotically flat exterior spacetime, for which there 
are fewer cases. The evolution within an asymptotically dS spacetime does not display 
any novel characteristics, and is discussed more briefly. The study of the ‘potential’ shows 
that there are two cases in which bubbles do not expand: either they start small enough, 
or their expansion is hidden behind a black-hole horizon (this possibility corresponds to 
A < 0). We will discuss in section 4.3 if any of these possibilities is realised in the simpler 
case of Higgs bubbles in an external flat spacetime, after the end of inflation. 

The standard evolution of sufficiently large bubbles is characterised by expansion, with 
their wall crossing the outer (dS) horizon. A typical example is presented in the Penrose 
diagram of £g. 7: the wall starts below the inner horizon and subsequently expands, pass¬ 
ing through both horizons and eventually reaching a speed close to that of light. The 
bubble grows in size and takes over part of the dS spacetime. The crucial question is 
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whether the bubble can engulf the total exterior spacetime, thus ending inflation. It is 
apparent from fig. 7 that this does not happen. Asymptotically the AdS bubble replaces 
only part of the spacelike surface r = oo, with the remaining dS space remaining unaf¬ 
fected. In other words, expanding bubbles are inflated away. Inflation precisely has the 
purpose of splitting the expanding universe into causally disjoint regions, and this limits 
the effect of the bubble growth: bubbles expand, but the dS space between them also 
grows. In the limit of inhnite inflation, both the bubbles and the exterior de Sitter phase 
acquire inhnite extent. 

We can estimate the asymptotic bubble size through the use of dS planar coordinates, 
commonly employed in the study of inhation. The metric has the form 

ds^ = —dtp + (^dVp + r'^dVt^ . (48) 

Assuming that the wall follows an almost null trajectory, we hnd that its location is given 
by 

Tp = rpo + ^ (l - . (49) 

Bubbles are created within the causally connected region, which extends up to 1/H at 
fp = 0. This means that its typical radius VpQ is of order IjH. Its subsequent growth 
extends this radius by 1/H}'^ It is reasonable then to expect that during inhation a typical 
bubble can be created with a certain probability within a causally connected region and 
then will roughly follow the general expansion of this region outside the horizon. It 
cannot, however, engulf the whole spacetime. The picture is completely diherent for an 
asymptotically hat exterior. As we shall see in the next section, in that case an expanding 
bubble can take over the whole spacetime. 

Another important question concerns the consequences for an outside observer of the 
AdS ‘crunch’ in the bubble interior. We discuss this issue in detail in appendix A.5. From 
the point of view of an observer deep inside the bubble, the coordinates in which the 
bubble appears as homogenous are those of an expanding and subsequently contracting 
open FRW universe with constant negative energy density. The bubble wall can be roughly 
identihed with the t = 0 surface in this slicing (see hg. 19 in the appendix). After a hnite 
(and short) time f, of order the AdS radius, a singularity forms in the bubble interior. 
However, this singularity never reaches the wall, as the latter expands with the speed of 
light. On the other hand, from the point of view of an external observer the bubble just 
expands forever (within either de Sitter or Minkowski spacetime). 

^^Note, however, that the physical bubble radius is obtained after multiplication by the divergent factor 
exp(i7tp). 
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Figure 8: Left; An example of how the dynamical evolution during pre-heating can bring 
back the Higgs field into the stable region before the instability takes over. The coloured 
lines show the potential at successive intervals of time, and the black line shows the tra¬ 
jectory of the Higgs field. Right; Maximal value of the Higgs field at the end of inflation 
(hend) that is brought back by a non-minimal gravitational coupling fu into the stable 
region, h < h^g^x, shown as a function of the Hubble constant during inflation and for 
different values of fn- 


4 Higgs evolution after inflation 

In this section we study the evolution of the Higgs field after infiation, considering that 
infiation ends with a matter-dominated phase, characterised by infiaton oscillations, fol¬ 
lowed by the reheating process, which ignites the usual thermal phase characterised by a 
gas of SM particles. 

4.1 Higgs evolution during pre-heating 

We start by considering the pre-heating phase, during which the energy density of the 
universe is dominated by the infiaton oscillations around its minimum. The interest of 
this phase lies in the case in which the Higgs potential has an extra mass term 
induced either by a non-minimal coupling to gravity fm? = or by a coupling to 

the infiaton 0 (m^ = In either case, the mass term rapidly shuts off after infiation. 

However, we will show that the induced m? can still have an important stabilising effect 
during the pre-heating phase. 
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Let us suppose that, at the end of inflation, the universe enters a matter-dominated 
phase, where the equation of state is that of a pressure-less gas, which is a good approx¬ 
imation when the inflaton held is oscillating before reheating. In this case the effective 
mass of the Higgs held is m? = and the Hubble rate scales as 

where we have set the value of the scale factor a at the end of inhation to unity. 

We consider a region in which, once inhation ends, the Higgs held has the value hend- 
If the Higgs mass term during inhation (rr? = is larger than (9/4)iL^ (he. 

if < —3/16), then the Higgs held is anchored at the origin and hend = 0. We are 
interested here in the opposite regime, when —3/16 < < 0 and quantum huctuations 

of the Higgs are generated during inhation; in this case, /lend is generally not zero. The 
subsequent evolution of the Higgs held is governed by the equation 


cPh 


, .dh 


+ 



d?h 5 dh a dV 
da^ 2a da H'^ dh 


Keeping only the mass term in V and neglecting the quartic term, we hnd {a/H‘^)dV/dh = 
—?i^Hh/a? and then the solution of eq. (50) is 


h\Cl) ^end Ci ^ 




k 


end ^ 




(51) 


where we have neglected the solution with h oc which is rapidly damped with 

respect to eq. (51), and where the last approximation is valid only for \^h\ ^ 3/16. 

As the amplitude of the Higgs held and the contribution to the potential from the 
term are both decreasing in time, we need to investigate if the Higgs held has time 
enough to roll down to the safe region h < /imax before the instability starts to become 
more important than the fading m? term, reverting the evolution of h. An example of the 
Higgs held behaviour is shown in hg. 8a. The time at which the instability starts driving 
the Higgs dynamics can be estimated by requiring that the quartic term in the scalar 
potential (A/i^/4 with A = — 61n/i^/Zi^g^i/e) is comparable with the mass term (m^/i^/2 
with at /^max, which corresponds to 


3 3 

(1^ = n ~- 

^max 7 7 2 

^'^max 


Thereby, the instability is avoided if 


(52) 


h 


end 


^ ^max*- 


max * 


(53) 


This estimate is conhrmed by the result of a numerical computation illustrated in hg. 8b, 
in which eq. (50) is solved using the full SM potential. 

recall that the Ricci scalar is given by it = —6(a/a -I- a? ja?) for a spatially flat universe. During 
inflation a = exp{Ht) with H constant, hence R = —12H^. In a matter-dominated phase (a oc we 
have R = —with = d/a oc while i? = 0 in the radiation-dominated phase (a cx H/^). 
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The qualitative conclusion, whenever 7^ 0, is the following. Higgs held values that, 


at the end of inhation, are even a factor of 0{2) above the instability scale h, 


max 


brought back into the metastability region {h < hmax) and saved from collapse into the 
AdS vacuum by the non-minimal gravitational coupling during pre-heating dynamics. 
The bound on the Hubble scale during inhation in eqs. (38) and (39) are correspondingly 
weakened by an 0(2) factor. On the other hand, in regions where the Higgs remains in 
the instability region > h max during pre-heating, the held quickly falls down into 

its deep minimum, in a time t ~ dTr/h max , unless a large enough temperature prevents 
the collapse, as we discuss in the following section. 

The dynamics discussed above is similar to the one in which the Higgs held is coupled 
to the inhaton held by a coupling of the form that generates a contribution 

to the Higgs mass. As discussed in section 3.1, if during inhation m > 3H/2 
Higgs inhationary perturbations are suppressed and the Higgs is efficiently anchored at 
h = 0. If instead m < 3H/2, Higgs huctuations are generated as in eq. (40) and they may 
pose a threat. After the end of inhation, when the inhaton held oscillates and its amplitude 
is redshifted away as 0 ~ the ehective mass squared of the Higgs, m? = 

decreases as a“^, exactly as the w? induced by the coupling {w? = One can 

therefore deduce the dynamics upon identifying the two m^. Of course, if the coupling of 
the inhaton held with the Higgs held is of a diherent nature, e.g. a non-renormalisable 
coupling of the form 0^/i^/Mpi, one needs to account for the diherent behaviour of the 
Higgs ehective mass. 

4.2 Higgs evolution during reheating 

In this section we study how the reheating process ahects the bounds on the Hubble 
constant H during inhation. Indeed, the dynamical evolution during the thermal phase 
can bring back the Higgs held towards the EW vacuum, even in regions where h has 
huctuated beyond the instability barrier {h > hmax) at the end of inhation. As a result, 
the bounds on H derived in section 3 are ehectively relaxed. 

At the end of inhation, the energy density of the universe is dominated by the coherent 
oscillations of the inhaton held 0 with energy density p^{t). The oscillations of 0, started 
at time t ~ l/TT, give a matter-dominated stage that gradually ends at t ~ l/T,^, where 
T^ is the inhaton decay width. The decay of the inhaton held into light degrees of 
freedom, which quickly thermalise via SM interactions giving rise to an energy density 
Pnit), initiates the radiation-dominated era of the universe. The process is described by 
the equations 



— —4:HrPR + T^p^, 



( 54 ) 
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Figure 9: Maximal value of the Higgs field at the end of inflation hend that is brought hack 
to the stable region, h < hma.x, o,s a function of the Hubble constant during inflation, for 
different values of and fn- This result, presented only in terms of ratios, negligibly 
depends on the absolute value of the instability scale /imax? which suffers from large uncer¬ 
tainties mainly due to the top guark mass. In the red region the inflationary fluctuations 
typically drive the Higgs to its negative-energy minimum, and therefore the corresponding 
values of hend and H are a highly unlikely outcome of inflation. In the parameter region 
below the dashed line, the field h rolls towards the SM vacuum, even in the absence of any 
thermal effect. 


where Hr = d/a = + pr) /{/iMf,/) is the time-dependent Hubble constant during 

reheating. 

The solution for the time evolution of p^ is 


pM 


P0(O) 

a3(f) 


P0(O) 


Stt 


(55) 


where the initial condition ^^(0) is given by the total energy density at the end of inflation. 
The second equation in the system (54) can be more conveniently written as 


dR _ 70^/2$ 

da + R/a ’ 


R = pro^ , 


<F = p^af , 


7 = 



rp2 

-^RH 


(56) 


where g^ is the number of degrees of freedom in the thermal bath (gf* = 106.75 in the SM) 
and Trh is the temperature of the system once all the inflaton energy is converted into 
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thermal energy at the decay time, 

( 57 ) 

Equation (56) can be approximately solved at the early stage of reheating (t <C E^^), 
by taking 1 in eq. (55) and neglecting the thermal-energy contribution to Hr 

{R/a <h). Once we express pr in terms of the effective temperature T, 

pad) ^ (58) 

the solution of eq. (56), at early times, gives the evolution of the temperature T (valid 
till the universe enters the radiation-dominated phase) 

T„.ax a-3/8(l - ' , (59) 

where kx = = 1.3 and k^ = = 0.54. 

The temperature T of the SM-particle gas raises from 0 to the maximum value T^ax as 
long as, soon after inflation, the scale factor of the universe a grows by an order-one factor 
in a time t ~ 1/H. After reaching Tmax, the temperature decreases as signalling 

the continuous release of entropy from the decay of the inflaton held. When this energy 
release ends, at time t ~ l/T^, the temperature is equal to Trh, which is called the 
reheating temperature, and then radiation cools in the standard way, T oc 1/a, due to 
space expansion. Note that the entire reheating process can be described by only two 
parameters, which we choose to be H and Trh- The decay of the Higgs condensate at the 
end of inhation has been discussed in ref. [26]. 

Let us now consider the evolution of the Higgs held throughout the thermal phase. 
Because of thermal corrections, the Higgs potential receives an extra mass term 
where rriT ~ gT and g represents the relevant combination of coupling constants. This 
expression for the thermal mass holds up to held values h < 2'kT. The thermal corrections 
to the potential can be approximated as [27] 

~ (o.21 - 0.0071 logio -^) , (60) 

where we added an exponential cut-oh at high values of h. Vr helps in stabilising the 
Higgs potential by shifting the instability to higher scales, in much the same way as the 
mass term due to the coupling does. 

Figure 9 shows the maximum allowed value of hend in order for the Higgs not to fall 
into its true vacuum at (or above) Planckian held values. A direct comparison with the 
right panel of hg. 8 shows that, for high enough reheating temperatures, thermal ehects 
are indeed of extreme relevance. 
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Figure 10: Minimal reheating temperature Trh needed to prevent the fall of the Higgs down 
into its deep true vacuum, assuming two different values for the instability scale /imax of 
the Higgs potential. 


It is not difficult to understand the behaviour of the maximum allowed value of hend 
as a function of the Hubble rate. Let us consider for simplicity the case f,H = ^ (left panel 
in £g. 9). From eq. (50) one can see that in a time scale of the order of the Hubble time, 
when the scale factor changes by order unity and the maximum temperature has been 
reached, the Higgs held changes by an amount (we neglect factors of order unity) 

(61) 

where we have approximated the zero temperature potential as Ah^/4 and one has to re¬ 
member we are considering the region where A < 0. Therefore, we obtain the approximate 
expression, valid soon after inhation, 

h - (^1 - ^ • ( 62 ) 

For h to roll towards the origin, a necessary condition is that the right-hand side of eq. (62) 
is smaller than one, which implies hend ^ Fmax/IA]^^^. This explains the approximate 
hatness of the curves in the right panel of fig. 9 for small H. For H 3> Tmax , the 
approximate scaling of the bound on hend as can be understood in the following 

way. Being the Hubble rate large, the term [a/H‘^)dV/dh in eq. (50) can be neglected 
up to the moment when the second time derivative term or the first time derivative term 
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become of the order of the potential term. This means that the Higgs field does not move 
much from its initial condition hend up to the moment when 


b dh a 2 (rp\i u rp 2 - 3 / 4 t 

This implies that the Higgs field starts moving away from hend when 

t2 \ 4/9 


a ^ = 


45 

8TL 


(63) 


(64) 


Imposing that at this value of a the hnite temperature term in the potential dominates 
over the negative quartic term gives 


h 


end 


< 

r\j 


/8^ 

, 1/6 ZJ-l/3rp4:/3 
\ max 

V45. 



(65) 


which reproduces the right scaling shown in £g. 9. 

In £g. 9 we let hend and H vary independently. However, the inflationary dynamics 
correlates the two variables, assigning a certain probability to hgnd for any given H. 
Although we have not used any relation between the two variables, in fig. 9 we have 
indicated in red the region in which the field h has overwhelming probability to slide 
towards large values and thus the corresponding parameters essentially cannot be the 
outcome of inflation. 

Figure 10 shows the minimal value of Trh for which the thermal corrections prevent 
the fall of h down to its deep minimum. In other words, it shows how the limit on the 
Hubble constant H corresponding to the orange area of £g. 6 can be relaxed, depending on 
the reheating temperature Trh- For sufficiently large Trh, regions in which the Higgs held 
huctuates around its instability scale hmax can be recovered by post-inhationary thermal 
ehects. On the other hand, the thermal phase cannot save regions in which the Higgs fell 
down into its deep negative-energy minimum during inhation. In the next section we will 
show that most of these regions eventually expand and hence a viable cosmology require 
that no such regions are produced during inhation. This excludes the red area in hg. 6. 


4.3 Bubble evolution in Minkowski spacetime 

The discussion of bubble evolution in an external Minkowski spacetime is analogous to 
the de Sitter discussion of section 3.3. A hrst diherence is that the ehective potential of 
eq. (46), which dictates the evolution of the bubble, is simplihed when we set iout = oo 
{i.e. (5 = 0). A second key diherence is that the external Minkowski space has no causal 
horizons: if, after inhation ends, bubbles expand at the speed of light, they engulf the 
whole space. 
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Figure 11: Penrose diagram describing an AdS bubble that expands in exter¬ 
nal Minkowski. Left: The wall trajectory corresponding to line C of fig. 15, in AdS 
space. Right: The wall trajectory corresponding to line C of fig. 15, in the Schwarzschild 
geometry. 

An important task is to determine whether Higgs bubbles expand or shrink. The 
complete analysis is presented in appendix A, where all the possible wall trajectories are 
determined. In summary, there are two scenarios in which bubbles do not take over the 
whole space: either they start small enough so that they shrink, or they expand but 
remain hidden behind a black-hole horizon (a possibility that corresponds to A < 0). 
In the following we examine if either of these possibilities is realised for Higgs bubbles, 
making them benign. 

We hrst consider bubbles with positive 

A = - {AnGay > 0, (66) 

which are bigger than their Schwarzschild radius and thereby can expand in the naive 
Newtonian way. This is the case depicted in the Penrose diagram in £g. 11, to be compared 
with £g. 7 for external dS. The black continuous curve denotes the trajectory of the wall: 
the bubble starts small and expands indehnitely within the asymptotically flat spacetime. 
The total space is constructed by patching the part of the diagram on the right of the 
wall with the part of the left diagram on the left of the wall. The shaded areas correspond 
to the parts that must be eliminated in order to join the remaining parts along the wall 
trajectory. From the point of view of an external observer, the bubble asymptotically 
expands at the speed of light and asymptotically reaches null infinity, filling all space. 
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Bubbles may also shrink because of their surface tension, if they are small enough 
{R < Rcr), and start with a small wall velocity. The critical radius, separating the two 
types of evolution, can be computed from the potential of eq. (46), but the resulting 
expression is not very illuminating. We present the complete discussion in appendix 
A.6. The result can be simplihed by assuming R = 0, and the Newtonian limit k -C 1 
(such that A > 0). In this case, the critical radius is obtained by extremising the sum 
of the surface and volume energy (dvri^^a — AnR^Vi^/S) with respect to R, thus hnding 
Rcr = 2a/Vin and a bubble mass M = IQna^/?)V-^. The exact result, valid even beyond 
the Newtonian approximation, is shown in £g. 12. In the ultra-relativistic limit the critical 
radius becomes Rcr = 3GM, slightly larger than the Schwarzschild radius 2GM. 

Figure 12 also shows the estimated M{R) corresponding to Higgs bubbles for different 
values of hin, the unknown Higgs value at its deep minimum. We use Idn ~ 
estimate the surface tension from the Newtonian expression 


cr 


dr 


1 / dh \^ 

2[^) +nh)-Vih,rr) 


hi 


-^ + Ar\X\hi> ^/\X\hl 


(67) 


minimised for a bubble thickness Ar 

The estimate for a is inserted in the full expression for M, eq. (90). The meaning of 
negative values of M is discussed in the appendix, in sections A.4 and A.6: they correspond 
to bubbles for which the negative volume contribution to the total energy budget is 
dominant. For our present purposes, the conclusion is that inflationary fluctuations create 
a number of bubbles with a variety of values of R and a, including bubbles with R > 
Rcr, which expand. Bubbles produced by inflationary fluctuations tend to appear with 
characteristic size i? ~ 1/if = iout and with negligible R. If the Higgs value is near the 
deep minimum of the potential, such bubbles are expected to have small or negative mass 
(because of the negative volume contribution to the energy) and, therefore, expand. 

Bubbles with sufficiently large surface tension a have A < 0 and their expansion is 
energetically disfavoured from a Newtonian point of view. However, there exist expanding 
solutions for such bubbles within general relativity. They are discussed in appendix A 
(see £g. 18). The crucial characteristic is that the expanding region is not accessible to an 
observer in the asymptotically flat space. In other words, an observer outside the bubble 
only sees a black-hole horizon, that protects him from its expansion. While such bubbles 
would be benign, A < 0 represents an extreme case: for the quartic Higgs potential, the 
surface tension gives a Planck-suppressed correction, such that 


A ~ 


M|i 


hi 


'^ir 


Mpi 


( 68 ) 


is negative only when the deep minimum of the SM potential is super-Planckian, hjn > Mpi. 
It is then impossible to make hrm predictions; strong gravitational effects may induce 
various dangerous effects. As long as the deep minimum is in the calculable sub-Planckian 
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region {hia<M-p\), bubbles have super-Planckian tension only for extreme Higgs field 
configurations, e.g. if the variation of h between the two minima happens within a sub- 
Planckian length. Inflationary fluctuations tend to create bubbles with bigger thickness 
and smaller surface tension. Furthermore, even if the condition A < 0 were initially 
satished, the bubble would evolve towards a smoother conhguration with smaller surface 
tension by reconhguring the Higgs held prohle. 

Our study has been carried out within 
the thin-wall limit, because this is the only 
setup for which an analytical treatment is 
possible. As we pointed out above, the re¬ 
alistic situation is more likely to involve 
conhgurations with a smooth transition re¬ 
gion from the interior AdS space to the 
false-vacuum exterior. For these, the fun¬ 
damental dynamics is mainly determined 
through the interplay between the nega¬ 
tive energy density in the interior and the 
positive contribution from the transition 
region. We expect that our analysis cap¬ 
tures the essential features of the evolution 
of such conhgurations as well. 



The important conclusion that we draw 
from our study is the following. No 
robust general-relativistic effect prevents 
large-field Higgs bubbles from expanding 
and engulfing all Minkowski space. As a 
result, a viable cosmology requires that no 
expanding bubbles are present in our past 
light-cone. In other words, the condition 
p{h —)• oo) < for A 60 e-folds must 


Figure 12: The boundary in the {R, M) 

plane that separates expanding bubbles from 
shrinking bubbles. The curves are the esti¬ 
mated masses of Higgs bubbles, as functions 
of their radius, for different values of Rn, the 
unknown Higgs field value at the deep mini¬ 
mum. 

hold and the red region in £g. 6 is excluded. 


4.4 The Higgs potential for 7 ^ 0 at zero temperature 

We conclude our study of the Higgs evolution after inhation with a remark concerning the 
non-minimal gravitational coupling As previously discussed, a coupling helps to 
stabilise the Higgs held during inhation. However, long after inhation, at zero temperature 
and H ^ 0, it could have an opposite ehect in the classical potential, in presence of a Higgs 
instability generated by SM interactions (Aeff < 0). In this setup an additional source of 
instability is generated by at Planckian values of the Higgs held, with no ehect on our 
discussion of the Higgs dynamics during the inhationary and post-inhationary phase. 
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SM potential Ve = AgeffC/if) h\lA for ^// < 0 



Figure 13: The SM Higgs potential Ve = y^Eea{hE)h%lin the presence of a negative fn 
coupling, written in the Einstein frame in terms of the canonical Higgs field He- 


In order to investigate the phenomenon, we focus on the real component h of the Higgs 
doublet = (0, h/\/2) and perform a Weyl rescaling to the Einstein frame gjf,^ = g^i, x / 
with / = /Mfy Then, the Einstein-Hilbert term becomes canonical and the action 

is 




\J det qe 




M, 


^'Re + zVH 
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/2 
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14(h) = 


14(h) 

V{h) 

P ■ 


+ 


(69) 

(70) 


We are studying the theory long after inflation, and therefore 17(h) = Aeff(h)h^/4/^ It is 
convenient to define a canonically normalised Einstein-frame Higgs held h^ through the 
equation dhE/dh = \fZ, where Z is given in eq. (70). The held h^; is such that He — h 
for h ^ Mpi and he —)■ cxd for h —)■ Mpi/ \/—^h (hence / —)■ 0). 

The Einstein-frame scalar potential becomes 


VE^hE) 


Aeff(h)h" 

4(l + errh2/M2;)2 


h=h(hE) 


(71) 


^^During inflation V contains an extra constant term Wj which dominates the energy density, and 
7/2 = is the Hubble constant during inflation. By expanding Ve at leading order in H/M^y 

for W / 0; we recover the Higgs mass term considered in the previous sections, m? = —= 
— The higher order terms were not relevant for our previous discussion. 
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In the limit of large He (which corresponds to h —>■ Mpi/ \/—^h), the denominator in 
eq. (71) nearly vanishes, while Aefr is negative. In that regime of held conhgurations, we 
hnd h{hE) — Mpi/ y/—^H oc Mpi exp{—^j2/3hE/Mpi) and the potential becomes 

VE{hE) ^ ( [S hE\ ^ ^ ^ 

forft^^oo. (72) 

Here Aefr is evalnated at h = Mpi/ yj—in) and is negative. The exponential behaviour 
in eq. (72) contributes to amplify the source of instability already present. The effect 
is shown in £g. 13, where we plot the effective coupling Agefr, dehned in analogy with 
previous effective quartic couplings by rewriting the potential in the Einstein frame as 
VE^hE) = AEeff(h£;)h|;/4. Negative values of appear to trigger a deep instability at 
Planckian held values h ~ Mpi/However, for nonzero the ultraviolet cutoh 
of the SM is no longer Mpi, but M-p\/\^h\ [28]. Hence, the instability just described 
takes place above that cutoh, where one is losing control of the theory [29]. It has been 
pointed out that small primordial black holes can seed Higgs-vacuum decay and enhance 
its rate [30]. However, this process crucially depends on the number density of black holes 
at a given epoch. 

5 The quantum gravity prediction for the Higgs mass 

In this section we explore how a speculative conjecture that has been put forward in 
the context of quantum mechanical completions of gravity can lead to a sharp correlated 
prediction for the Higgs and top-quark masses. The intriguing result is that this prediction 
agrees quite well with the measured values of these masses. The reasoning is essentially 
based on two points. 

1. The empirical observation that we live in an accelerating universe. 

2. The theoretical conjecture that quantum gravity is ill-dehned in de Sitter space [31, 
32]. 

The difficulties with dS quantum gravity have been argued from various perspectives [31, 
32] . We summarise here in a very schematic way some of the arguments against a stable 
dS space, reviewed in [33], extending them in light of some more recent developments. 
There is no positive conserved energy in dS (and, as a consequence, there cannot be unbro¬ 
ken supersymmetry). There is no classical compactihcation of ten- or eleven-dimensional 
supergravity to dS space, and stable dS space cannot be obtained from any string or 
M-theory. Even in the general setting of quantum gravity, beyond the particular UV 
completion ohered by string theory, other problems arise. It has been suggested that the 
quantum Hilbert space in dS is of hnite dimension, limiting the variations of complex 
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constructions. Given that the Gibbons-Hawking temperature sets a minimum temper¬ 
ature, the hnite dimensionality of the Hilbert space sets a maximum time scale, the 
so-called recurrence time [32], In particular, this leads to the problem of the so-called 
Boltzmann brains [34], and it has been suggested that its resolution calls for an unstable 
universe [34], More generally, a rigorous dehnition of the Hilbert space in dS seems to be 
problematic [31]. In quantum gravity, it is difficult to dehne precisely local observables 
and one can rely only on asymptotic quantities, such as the S-matrix in Minkowski space 
and the boundary correlators in AdS. However, in dS, where asymptotic states fall behind 
the horizon, no such precisely dehned observables seem to be present. 

In addition to these problems that have been known for some time, it has been found 
in [35-37] that there is a sharp universal bound on how much reheating volume slow-roll 
inflation is capable to create without being eternal. This is given by where S^s 

is the entropy of the would be de Sitter space with Hubble rate evaluated at the time 
of reheating. Larger overall expansion is possible only by making space inhnite. This 
generalises at the quantum level the bound on the duration of inflation found in [38], 
and determines the universality of such a bound under the number of helds involved, 
higher derivative corrections, number of space time dimensions and slow roll parameters. 
The same phenomenon, i.e. not being able to produce arbitrarily large hnite volumes, is 
shared by false vacuum inhation. These results seem to suggest that there are bounds on 
the kind of global spacetime structures that local quantum held theory can generate. This 
is related to, and somewhat supporting, the argument against de Sitter space in nature 
that we follow here. 

Though none of these arguments raise to the level of a proof of the inconsistency of dS 
space, they clearly give an idea of the conceptual difficulties that arise when considering 
quantum mechanics in dS space. 

It could well be that the problems of asymptotic dS space are circumvented by Planck- 
ian dynamics, which can for example open channels for vacuum tunnelling to a ‘landscape’ 
of other minima with zero or negative vacuum energy. This is certainly a possibility. 
However, it is interesting to note that, even without any special hypothesis about the 
gravitational sector, the SM Higgs offers an easy way out to the problem. A solution is 
automatically found if the present dS space is only metastable. 

As soon as the decay rate per unit space time volume T of the false vacuum is non 
zero, the asymptotic space is not dS. True vacuum bubbles are nucleated and expand at 
the speed of light. There are two critical values of the decay rate [39], both valued around 
with H\ being the Hubble rate of the would-be de Sitter region. 

If T is larger than the largest critical point, bubbles percolate, £11 the whole space, 
inflation ends globally, and the asymptotic spacetime is the one of the true vacuum (if 
this is AdS, the instability grows and leads to a singularity in about one Hubble time). 
The value of this critical point is T 2 /H^ = 9/47r ~ 0.71 [33]. 


37 


The precise value of the second critical point is unknown, but is bounded to be in the 
range [39] 1.1 x 10“® < Ti/H\ < 0.24. If the decay rate lies between these two critical 
points, then bubbles percolate and long chains of bubbles form, connecting arbitrarily 
distant point in the otherwise dS space. Finally, when the decay rate is even slower than 
the second critical point, then bubbles do not percolate, and the asymptotic spacetime is 
the one called ‘false vacuum eternal inflation’. In this phase, there is an inhnite amount of 
space that keeps inflating, but each single point decays at some time into the true vacuum 
through the nucleation of a bubble.^® 


On the other hand, there is no question that AdS and Minkowski are well-dehned 
spaces from the point of view of quantum gravity, although we have a non-perturbative 
formulation of quantum gravity only for AdS. So, it is possible that quantum gravity 
allows for an eternally inflating false-vacuum space-time, but not for dS space. We will 
show that the opportunity of circumventing dS space through the Higgs held is given to 
nature only for a narrow range of Higgs boson masses. Interestingly, it seems that nature 
did not miss the opportunity because the measured Higgs boson mass lies exactly within 
this range. 

Let us for a moment accept the quantum-gravity arguments against stable dS space, 
and let us assume that our universe, which we observe to be today in a dS phase, escapes 
the problems through a future decay of the Higgs vacuum. This implies the bound on the 
Higgs boson mass [1] 


Mh 

GeV 


< 129.6 + 2.0 



173.34 


0.5 


a3(Mz)-0.1184 \ 
0.0007 ) 


(73) 


The vacuum decay rate induced by the SM Higgs instability, exponentially suppressed 
by the action of its bounce solution, S ~ 87r^/(3|A|) ~ 10^, is always faster, and typically 
much faster, than what is needed to avoid the dS problems. Indeed, for comparison, the 
de Sitter entropy is S'ds 10^^°. In general, the action of any Coleman-de 

Luccia instanton out of a false vacuum with positive energy density is never larger than 
S'ds, no matter how high we make the false vacuum barrier (see the discussion of [33] 
for a review).^® Therefore, when the Coleman-de Luccia instanton is present and can 
be reliably computed, the lifetime is bound by the Poincare recurrence time of de Sitter 
space, of order up to logarithmic factors. It is expected that this lower bound on 

the decay rate is a property shared by all theories with a de Sitter false vacuum within 
held theory and perturbative gravity. 

^^Bubbles continuously form and collide an infinite number of times among each other but they do so in 
relatively smaller and smaller regions: there are points in the inflationary space for which the probability 
that they are connected by a stream of bubbles is zero. When bubbles collide in our past light cone, a 
very sharply defined disk shape is impressed in the CMB [40]. The optimal analysis to search for such a 
signal in the WMAP data has been recently performed, with no evidence found [41]. 

^®One can check this by taking the limit in which the wall energy density goes to infinity in eq. (3.16) 
of the paper by Coleman and de Luccia [7]. 
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There are also lower bounds on the Higgs mass. A hrst bound is obtained by requir¬ 
ing that the tunnelling rate away from the EW breaking state with small and positive 
cosmological constant is not faster than the age of the universe. In the current universe, 
the volume of our past light cone at the current time Tu is 


Vol 4 (r„) = OMHl , 


(74) 


where = A^/3Mp^ is the Hubble rate produced by the observed vacuum energy A^. 
Imposing that the probability p = e“^°* 4 (T'c/)r universe experienced vacuum decay 

in the past is small enough, implies an upper bound on the vacuum decay density rate T 


^ Voi4(r[,) ip 


(75) 


and, within the SM, a lower bound on the Higgs mass [1] 
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A stronger lower bound is obtained from the requirement that the universe underwent 
a hot phase. There are good reasons to believe that the universe has been very hot 
at an early epoch. Indeed, processes such as inflation and leptogenesis suggest that 
the primordial universe reached high temperatures. We have seen in section 4.2 how 
a large reheating temperature helps in forcing the Higgs to its weak scale meta-stable 
minimum. However, such high temperatures could have prematurely destabilised the 
Higgs metastable vacuum. The requirement that this did not happen implies 
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_______ 1 


±1. (77) 


Equations (73) and (77) dehne a fairly narrow range of possible Higgs masses (see £g. 14). 
Loosely speaking, one might claim that quantum-gravity favours 


122 GeV <Mh< 129.4 GeV for Mt = 173.34 GeV. (78) 


Given that the Higgs mass is now precisely measured, one can better use Mh as input and 
predict the top mass in the range 

(79) 

The coincidence that the Higgs and top masses are within the predicted range can be 
viewed as an indirect indication that nature took the opportunity offered by the Higgs to 
avoid the problem of an asymptotic dS space. 


171 GeV <Mt< 175 GeV. 
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Figure 14: The allowed meta-stahility window of the Higgs mass. The ellipse indicates the 
measured values of Mh and Mf. The orange region is excluded by assuming that “stable 
dS” is unacceptable. The red region is excluded by vacuum decay at zero temperature. 
The pale-red region is excluded by the requirement that the universe must have been hot 
in the past (the dashed red curves show boundaries for different values of the reheating 
temperature). The bottom panel shows the same result in the full range of a-priori possible 
Higgs masses, in order to emphasise the smallness of the surviving meta-stability region. 


6 Conclusions 

Assuming that the SM holds up to large energies, we studied under which conditions the 
cosmological evolution does not disrupt the electroweak vacuum, in spite of the presence 
of an instability of the SM effective Higgs potential V{h) at held values h > hmax- 

As a preliminary step, in section 2 we clarihed the gauge-dependence of the effective 
potential. The Nielsen identities show that the gauge-dependence of the effective action 
corresponds to different ways of parameterising the same physics in held space: the phys¬ 
ical content of the ehective action is gauge-independent. We have shown how this implies 
that the classical equation of motion (as well as the related Langevin and Fokker-Planck 
equations used later) are gauge-independent, because the gauge-dependence of the ehec¬ 
tive potential is compensated by the gauge-dependence of the kinetic term. Furthermore, 
we showed how, in the basis in which the kinetic term is canonical, the full ehective po¬ 
tential becomes gauge-independent in the limit in which only the leading-log corrections 
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are retained (which, for our purposes, is a very good approximation, see fig. 1). For the 
present best-fit values of the SM parameters one has hmax ~ 5 x 10^° GeV, but hmax can 
vary by orders of magnitude if the top mass Mt is varied within its uncertainty band. 

Next, in section 3 we studied Higgs fluctuations during inflation. In our study, we 
also took into account the effect of a Higgs mass m? induced, during inflation, either by 
a mixed quartic coupling between the Higgs and the inflaton, or by a non-minimal Higgs 
coupling to gravity ^h- Not being radiatively stable, such a coupling is expected to be 
generally present and leads to m? = where H is the Hubble constant during 

inflation, given by FT ~ 8 x 10^^ GeV-y/r/0.1. Present cosmological data constrain r < 0.1, 
but future measurements will have greater sensitivity. 

If < 9iF^/4 the Higgs undergoes inflationary quantum fluctuations, which we com¬ 
puted via a Langevin equation that bypasses the need of imposing appropriate boundary 
conditions encountered in the Fokker-Planck equation used in previous works. We find 
that the parameter space in the plane iF/hmax vs splits into 3 regions (see fig. 6): 

• ‘green’ region, where the Higgs remains below its instability scale at the end of infla¬ 
tion, and thus inflationary fluctuations do not destabilise the electroweak vacuum. 

• ‘orange’ region, where the Higgs can probe field values above the instability scale 
(|F| > Fmax), but quantum fluctuations dominate over classical evolution and pre¬ 
vent the Higgs from falling into its true AdS minimum; the ultimate fate of the 
Higgs is determined by post-inflationary dynamics. 

• ‘red’ region, where the Higgs fluctuates above the instability scale and falls down 
into its true minimum, presumably ending inflation in that patch of space. 

In section 4 we followed the evolution of the Higgs field through the reheating process, in 
order to assess the viability of parameters corresponding to the ‘orange’ region. Thermal 
effects can rescue the Higgs field, letting it slide towards the origin of the SM potential, 
if the reheating temperature after inflation Trh is sufficiently large. We derived upper 
bounds on Ff/Fmax, for given Trh, as shown in fig. 10. The result is that thermal effects 
can easily make the ‘orange’ region cosmologically acceptable. 

On the contrary, we found that the ‘red’ region is problematic. By approximating 
the large-held Higgs patches as spherical bubbles with small thickness, we could perform 
a general relativistic computation in order to determine whether such bubbles shrink or 
expand. The computation addresses several relevant counter-intuitive phenomena. While 
we identihed mechanisms that can make some of the bubbles innocuous (small bubbles 
with low wall velocity shrink, bubbles with large tension expand hidden behind a black- 
hole horizon), we hnd that inhation produces Higgs ‘bubbles’ that expand, at least as long 
as they are in the computable sub-Planckian regime. During inhation these bubbles are 
not lethal, as they remain behind a de Sitter horizon and are diluted by space expansion. 
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However, after inflation they keep on growing at the speed of light, eventually swallowing 
all space. Therefore, we must require that inflationary fluctuations do not produce any 
of these regions in our past light-cone. 

This leads us to our hnal result: the ‘red’ region of £g. 6 is excluded. If \^h\ < 0.01 one 
needs a Hubble constant smaller than 0.045 hmax- This constraint gets weaker (stronger) 
for negative (positive) ^h- e-d- H < 10^ hmax for ~ —0.03. A small negative 
however leads to a new, super-Planckian instability of the SM potential in the Einstein 
frame, see £g. 13. In a similar way, a direct coupling of the inflaton to the Higgs could also 
relax the limits on H but, contrary to the case of h does not lead to any instabilities 
at large held values. 

Finally, in section 5 we explore a new speculative idea. Assuming that the present 
acceleration of the universe is due to a small cosmological constant, and accepting the 
conjecture that quantum gravity is ill-dehned in a de Sitter space, we argue that vacuum 
decay is a necessary way out for the universe. We show that vacuum decay triggered by 
the Higgs instability is fast enough to resolve this conceptual problem. 

Basically the SM phase diagram in the (M^, Mh) plane is reinterpreted: the instability 
region remains ‘bad’, the stability region becomes ‘bad’, and the only ‘good’ region is the 
narrow meta-stability strip of parameter space. As discussed in section 4.2 a large enough 
reheating temperature may play an important role in the universe, and the requirement 
that thermal effects do not induce an excessively fast vacuum decay provides a further 
restriction in the Higgs and top masses, as shown in fig. 14. One could view this restriction 
as a remarkably precise post-diction for the Higgs or top masses. 
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A Evolution of bubbles 


In this appendix we study the evolution of a region of true vacuum with negative vacuum 
energy density, which lies within the false-vacuum asymptotically flat or de Sitter space. The 
basic question is whether this region, which we call a bubble (assuming spherical symmetry), 
expands or contracts. For an outside observer, the presence of the bubble has a gravitational 
effect equivalent to the presence of a central mass. As a result, the exterior metric is of the 
Schwarzschild or Schwarzschild-de Sitter (SdS) type. We study an idealised configuration with 
constant vacuum energy density in the interior and exterior of the bubble, as well as constant 
surface tension. The study of a realistic bubble, corresponding to a space-time dependent Higgs 
configuration, is not possible analytically. However, we believe that our treatment captures the 
main aspects of the problem, determined essentially by the difference in the local energy density 
of the Higgs field on either side of the bubble wall. 

We employ the thin-wall approximation and parameterise the wall and the inner and outer 
space as described in section 3.3. The interior of the bubble is assumed to be a part of anti-de 
Sitter spacetime, described by the metric (41). The exterior of the bubble is described by the 
metric (42). The case Hout 7^ 0 corresponds to the SdS spacetime, while the case V),ut = 0 to an 
exterior Schwarzschild metric. The metric on the wall is given by eq. (43). 

A.l Matching the geometries 

The metric must be continuous over the whole space. This means that 


fin{R)V = 

ei [r^ + fin{R)y^\ 

(80) 

foutiR)i = 

62 (R^ + font{R)y^\ 

(81) 


where ei = ± 1 , 62 = ±1 are possible sign choices and a dot denotes a derivative with respect to 
r. Since /in > 1, the value of ei determines the relative flow of the two timelike coordinates rj and 
r. It is natural to make the choice ei = 1, which is also the only consistent choice (see below). 
We consider only this value in the following. The relation between t and r is more complicated 
because font can be negative. We follow the convention of [42], according to which the flow of 
proper time is such that future-directed world lines correspond to a growing Kruskal-Szekeres 
coordinate V (so that H > 0). 

The matching of the two regions can be done following [42]. The four-velocity of a point 
on the wall is = (?), R, 0) and = (i, R, 0) in each of the frames. A spacelike vector 
perpendicular to the wall must be orthogonal to U^. In order to determine it uniquely, we 
have to specify whether it points towards the interior or the exterior. For spaces with horizons, 
such as the exterior space, we adopt the convention of [42]. We assume that points towards 
increasing values of the Kruskal-Szekeres coordinate U. We also assume that the exterior lies 
on the ‘right’ of the wall in the Penrose diagram. For the AdS space in the interior, we assume 
that points towards increasing values of the global coordinate r. With these conventions 
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points from the interior towards the exterior. It is given by 


Ct = f /out i, o) = f e2(/out + R^Y'\ o) , (83) 

\/out J \Jout j 

in each of the frames. It has been normalized to —1. 

The junction conditions connect the discontinuity in the extrinsic curvature to the surface 
tension: 

(Xout))- - = -47r(TG<5V (84) 

We match the QQ component of the extrinsic curvature (the other components give equivalent 
relations), which is Kqq = Evaluated on either side of the wall, it is given by 

{Kin)e0 = (/in + = Ani?, (85) 

{Kont)ee = e 2 {font + R^f^^R^PoutR. ( 86 ) 

Thereby, the 69 matching condition is 

e2(/out + ^')^/' - (/in + R^ f^^ = Pout - Pin = -^irGaR. (87) 


A.2 Bubbles in asymptotically flat spacetime 


We consider first the case lAut = 0, which corresponds to an exterior Schwarzschild metric. We 
shall discuss later the case Eout A 0, corresponding to a SdS spacetime. 

The square of eq. (87) can be put in the form: 


2GM 




R^ + 2e2KR^ 



2GM 

R 



( 88 ) 


For large R and non-relativistic wall velocity, the last parenthesis becomes equal to 1. The 
resulting expression indicates that the mass M of a large-radius bubble is dominated by a 
volume contribution proportional to k? — The total volume effect can be negative or 

positive, depending on the value of 

e = sign . (89) 

As a result, it is possible for the total mass M to become negative. 

Solving eq. (87) for M one finds a result that, for €2 = 1, has a simple Newtonian interpre¬ 
tation: 

/ r2 . \V2 

R^ -|- 2kRP ( ^ £2 —^ R^ j 1 ( 96 ) 

with K = AttGct. For small i?, the mass M attributed to the bubble of AdS by an outside observer 
contains a volume term proportional to — The contribution —1/Pf^ corresponds to the 


2GM = - 


il 
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vacuum energy density, while —k? reproduces correctly the gravitational self-energy of the wall. 
The second term in eq. (90) can be expanded for small R and R. One recovers the surface energy 
of the bubble, with nonrelativistic correction, and the surface-volume binding energy [39]. The 
leading term for small R is the positive surface energy ~ kR? , which indicates that small bubbles 
tend to collapse in order to minimise their energy. The case £2 = — 1 does not lead to solutions 
with a simple Newtonian interpretation, even though it contains acceptable configurations for 
the global geometry. 

By squaring eq. (87) a second time, we can express the ‘kinetic energy’ R? in terms of a 
conserved ‘energy’ E and an effective ‘potential energy’. We express the result as the equation 
for the one-dimensional motion of a particle in a ‘potential’ V 


dR 

df 


+ V{R) = E, 


where 


and 


V{R) = - 


'em + eR^\ 7 ^ 

I G2MV 




em = sign(M). 

The dimensionless ‘coordinate’ variable R and the ‘time’ variable f are defined as 


2k 

R = pR, f = -^T. 

T 


The parameter p, defined as 


1 


P = 


1 


— K 


(91) 

(92) 

(93) 

(94) 


2G|M| 

sets a characteristic inverse length-scale, while 7 parameterises the surface-energy term in V: 

2k 1 -'2 


(95) 


7 = 




1/2 ■ 


2 1 7 

i.e. K = 


472 + 46 - 


(96) 


The form of the solutions of eq. (91) can be revealed more easily through the following 
observations: 


• The sign 62 disappeared when performing the second squaring, so that eq. (91) describes 
the solutions of eq. (87) with both values of 62 . We can rewrite eq. (87) in terms of the 
new parameters as 

An = /3out + 47rGaR = ^ + eR^ + ^R^'^ , (97) 

where we have used eq. (90). For positive-mass bubbles (e^ = 1) we have An > 0. This is 
obvious for 6 = 1. It also holds for e = — 1 , because 72 > 4 in this case. We conclude that 
the only consistent value for ei for positive-mass bubbles is ei = 1 (the value we assumed). 
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• We can also write 



(98) 


from which it is apparent that, for positive-mass bubbles, / 3 out is positive and 62 = 1 for 
e = 1 , while /3out changes sign at -R = 1 for e = — 1 . 


• For negative-mass bubbles (M < 0 or Cm = —1) the variable t is always timelike. The 
value of t 2 determines the relative flow of t and r. It is natural to make the choice 62 = 1 


in this case. The possibility €2 = — 1 does not lead to a physical solution, as we discuss in 
subsection A.4. 

• It is apparent from eqs. (92), (95) that, for fixed and k, the total energy S is a function 
of M. As a result, the nature of the various solutions of eq. (91) is directly related to the 
mass of the bubble. 

• The ‘potential’ is maximal at i? = i?max) given by 



(99) 


The value of the ‘potential’ at its maximum is 



F(i?max) = -3 


( 100 ) 


'max 


For positive-mass bubbles (e^ = 1) we have i?max > 1- 

• The Schwarzschild radius of a bubble with positive mass M is rjj = 2GM, which, in terms 
of the variable R, becomes 



( 101 ) 


Rh 


This relation determines the location of the horizon on a solution of eq. (91) with given 
E. Making use of the definition (92) of the ‘potential’, we can write 



( 102 ) 


For e = —1 the curve — 7 ^/i?//, depicting the location of the horizon, is tangent to the 
curve V = V{R) at iZ = 1. For e = 1 the curve for the horizon is always located above 


the curve for the ‘potential’. 

The above features are depicted graphically in fig. 15. The solid black curve depicts the 


‘potential’ V{R), which has a maximum at R = i^max- The dashed blue curve indicates the 
location of the horizon. For e = 1 (i.e. 1/-^^^ > the curve for the horizon is always located 

above the ‘potential’. For e = —1 (i.e. l/^f^ < the curve for the horizon is tangent to the 

‘potential’ at R = 1. The function /3out changes sign at this point. In the centre plot we have 
separated with a red vertical dashed line the regions in which e (and /3out) has opposite signs. 


max* 
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Figure 15: The ‘potential’ of eq. (92) for 7 = 8 and e = 1, = 1 (left), for e = —1, 

em = 1 (middle), and for e = 1 , = —1 (right). 


The various types of trajectories can be deduced from these plots. We plot a few lines with 
constant E that stop when E = V: at this point R = 0 and the motion of the wall is reversed. 
There are various types of trajectories for which the bubble expands indefinitely. If e = — 1, 
such evolution can be obtained only for 62 = — 1 . 

A.3 Evolution of positive-mass bubbles 

We consider first the case M > 0, or, equivalently, Cm = 1- The evolution of the wall is 
best depicted using Penrose diagrams. The diagrams for the most characteristic types of wall 
evolution are presented in hgs. 16, 11, 18. Each hgure contains a pair of diagrams. In each 
pair, the left diagram depicts AdS space, which has the simple structure of a cylinder, with the 
line marked r = 0 corresponding to its centre and the line r = 00 to conformal infinity. The 
right diagram represents the complete Schwarzschild geometry, which includes two singularities, 
marked r = 0, and the corresponding horizons. The black continuous curve in each diagram 
denotes the trajectory of the wall. Thick black lines denote singularities, dashed lines horizons 
and dotted lines denote conformal infinities. The total space is constructed by patching the part 
of the left diagram on the left of the wall with the part of the right diagram on the right of the 
wall. The shaded areas correspond to the parts that must be eliminated in order to join the 
remaining parts along the wall trajectory. 

The crucial relation for the fate of space is between the gravitational self-energy of the wall 

and the vacuum energy Naively, one expects that, if (i-®- ^ = !)> large 

bubbles will grow indefinitely because the system gains energy in the process. In the opposite 
case with < nf (i.e. e = —1), the bubbles will shrink for similar energetic reasons. These 
simple expectations, which are based on Newtonian intuition, are only partly fulfilled in the 
complete analysis. More complicated scenarios are realised as well. 
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Figure 16: Small bubbles with small initial wall velocity do not expand. Left: 
The wall trajectory corresponding to line A of fig. 15, in AdS space. Right: The wall 
trajectory corresponding to line A of fig. 15, in the Schwarzschild geometry. 

Case e = 1 

Small babbles with small initial wall velocity do not expand 

Line A of fig. 15 describes the evolution of a bubble whose volume energy receives its largest 
contribution from the negative vacuum energy density > nf). However, the surface 

contribution to the energy, arising from the wall tension, is the dominant factor and tends to 
make the bubble shrink. The bubble has small initial wall velocity, which prevents it from 
evolving to a size sufficiently large for the volume contribution to the energy to dominate. As a 
result the surface tension wins: the bubble reaches a maximum size and subsequently collapses 
falling within its own horizon. 

The space corresponding to this solution is depicted in fig. 16. It results from eliminating 
the shaded areas in each of the two Penrose diagrams and patching the remaining parts along 
the wall trajectory. 

Small bubbles with large initial wall velocity expand 

Line C of fig. 15 corresponds to a bubble with similar characteristics as in the previous case, but 
with much larger wall ‘kinetic energy’. This is is apparent by the fact that the total energy is 
less negative. We can consider a bubble that starts very small (with almost vanishing r ox R). 
Even though the surface contribution to the ‘potential energy’ dominates, the initial velocity is 
sufficiently large for the bubble to expand. Eventually the bubble develops a size for which the 
volume contribution to the ‘potential energy’ becomes dominant over the surface contribution. 
From this point on, the bubble expands indefinitely, with its wall approaching asymptotically 
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Figure 17: The evolution of negative-mass babbles. Left: The wall trajectory cor¬ 
responding to line H of fig. 15, in AdS space. Right: The wall trajectory corresponding to 
line H of fig. 15, in the negative-mass Schwarzschild geometry. 

the speed of light. 

The corresponding evolution of space is depicted in hg. 11. After a hnite time g the wall 
reaches the boundary of AdS space. As the AdS boundary is timelike, there is a Cauchy horizon, 
beyond which the spacetime cannot be determined without additional boundary conditions. It is 
expected that, within the full theory beyond the thin-wall limit, a spacelike singularity develops 
before the Cauchy horizon [24]. This is depicted by a thick blue line in the left diagram of hg. 11. 
From a mathematical point of view, the solution also describes the reverse process. 

Large bubbles expand 

Line B of hg. 15 describes the evolution of a bubble so large that its surface tension is irrelevant. 
The bubble starts with inhnite radius, shrinks to hnite size and then re-expands. There are 
two singularities in the Penrose diagram of AdS space, starting from the points at which the 
wall trajectory reaches the boundary [24]. The whole trajectory lies with the region I of the 
Schwarzschild geometry. 

Case e = — 1 

Small bubbles with small initial wall velocity do not expand 

Line D of hg. 15 describes evolution very similar to that for line A. The contribution from the 
surface tension dominates the ‘potential energy’, while the ‘kinetic energy’ is small. The bubble 
expands up to a certain size, and subsequently recollapses. The space is described by Penrose 
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Figure 18: Large bubbles with expand behind the horizon. Left: The 

wall trajectory corresponding to Line G of fig. 15, in AdS space. Right: The wall trajectory 
corresponding to Line G of fig. 15, in the Schwarzschild geometry. 

diagrams very similar to those of fig. 16. 

Line E of fig. 15 describes a similar scenario, but now the extrinsic curvature fiout (or, 
equivalently € 2 ) changes sign during the evolution. This implies that the wall trajectory crosses 
regions IV, III and II of the Schwarzschild geometry instead of the regions IV, I, II (see fig. 16). 

Small bubbles with large initial wall velocity expand behind the horizon 

As we have seen already, the case e = — 1 may lead to evolution that cannot be deduced through a 
purely Newtonian approach. For < nf the wall self-energy dominates the negative vacuum 
energy, so that the growth of the bubble seems energetically unfavourable. The Newtonian 
intuition suggests that such bubbles cannot expand. However, there is a relativistic solution 
described by line G of fig. 15. The corresponding space evolution is depicted in fig. 18. The 
crucial difference with respect to the case e = 1, depicted in fig. 11, is that the wall trajectory 
is located within the regions IV and III, instead of the regions IV and I of the Schwarzschild 
space-time: in simpler words the bubble expands inside its Schwarzschild radius. The extrinsic 
curvature fioutR of eq. (86) changes sign along the trajectory G of the wall, while it stays positive 
for the trajectory C, as can be seen in fig. 15 [42]. Asymptotic regions of flat space-time survive: 
the growth of the AdS region and its singularity are hidden behind the horizon and do not affect 
an observer located in region I.^^ 

^^A pictorial representation of an analogous situation for a dS bubble is given in fig. 13 of [42], in which 
case the AdS singularity is absent. 
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Large bubbles expand behind the horizon 

Line F of fig. 15 describes a large bubble that initially shrinks, reaches a minimal size and 
subsequently expands. The whole evolution lies entirely within the region III of the Schwarzschild 
space-time and is hidden behind a horizon for an observer located in region I. 

A.4 Evolution of negative-mass bubbles 

We next turn to the solutions with negative mass M, or, equivalently, Cm = —1- The metric 
(42), with Wut = 0, has a naked timelike singularity at r = 0 in this case. However, this metric 
is relevant only for the bubble exterior, while the interior is described by the AdS metric (41). 
As long as the bubble expands and the wall moves to increasing values of R, the global geometry 
is free of singularities. 

The form of the ‘potential’ for 1/1’^ > M < 0, depicted in the right plot of fig. 15, allows 
for such a solution. For e = 1, = — 1, the ‘potential’ has a positive maximal value. On the 

other hand, the ‘energy’ E of eq. (92) is always negative. This allows for only one possible type 
of solutions, the one corresponding to line H of fig. 15. It represents a bubble that starts with 
infinite radius, shrinks to a finite value of R, and subsequently re-expands. Its mass is negative, 
because the radius is always sufficiently large for the negative volume contribution to the energy 
content to be dominant. The Penrose diagram for this solution is depicted in fig. 17. There are 
two singularities in the Penrose diagram of AdS space, starting from the points at which the 
wall trajectory reaches the boundary. The Schwarzschild metric with negative mass has a naked 
singularity at r = 0, depicted by the vertical solid line in the right plot of fig. 17. However, this 
singularity is irrelevant for our problem because it is eliminated when the white areas of the two 
plots are joined along the wall trajectory. 

It must be pointed out that it is not possible to construct negative-mass solutions corre¬ 
sponding to horizontal lines extending from R = 0 to the ‘potential’ in fig. 15. As can be seen 
from eq. (98) such lines would require /3out < 0 , or, equivalently 62 = —1. As we have already 
remarked, this choice would require the timelike coordinates r and t to increase in opposite 
directions. For positive mass, the Schwarzschild geometry has sufficient structure to permit 
solutions with both signs of /3out) such as the one corresponding to line G of fig. 15, which is 
depicted in fig. 18. However, for negative mass, the Penrose diagram cannot be extended beyond 
that depicted in fig. 17, unless completely disjointed regions are introduced. For this reason, the 
only meaningful solution is the one of fig. 17. 

A.5 The AdS ‘crunch’ 

As we have seen, the particular structure of AdS implies that the evolution of the bubble must 
lead to a singularity. This is apparent in fig. 11: the form of spacetime, after the finite time 
r] at which the wall reaches the timelike boundary, cannot be determined without additional 
boundary conditions. It is expected that, for a physical system that realises an approximation 
of the idealised bubble evolution that we consider, a spacelike singularity must develop in the 
interior of the bubble [24]. 
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For the problem at hand, the nature of this singu¬ 
larity can be understood through the picture of the AdS 
‘crunch’ presented in [7]. As shown there, a part of AdS 
space can be viewed as an open Friedmann-Robertson- 
Walker (FRW) universe with negative energy density. 

The coordinate change [43] 

r = ^in sin — sinh ?/) 

^'m 

i ( r‘^\ r? 

cos— = U + ;^j cos— (103) 

puts the AdS metric (41) in the form 

ds^ = —dP -b sin^ (dip'^ + sinh^ V’ dO^') ■ (104) 

-tin ^ ' 

This metric describes an homogeneous FRW universe 
that is born with a big ‘bang’ at t = 0 and collapses 
in a big ‘crunch’ at t The coordinates t, V’ do 

not cover the whole AdS space, but only a triangular 
patch of its Penrose diagram. This is bounded by the 
dot-dashed null lines in fig. 19. Even though these lines 
represent only a coordinate singularity in the idealized 
picture of a pure AdS bubble interior, the big ‘crunch’ 
becomes a true physical singularity in the presence of a 
fluctuating Higgs field, as argued in [7]. From continuity 
it is apparent that the Higgs fluctuations become very large in the neighbourhood of the null 
line. 

The connection of the AdS ‘crunch’ to the bubble evolution can be obtained by establishing 
the relative position of the FRW ‘triangle’ and the wall trajectory on the Penrose diagram. As 
the FRW observer views a homogeneous universe, the AdS patch to which he has access must 
be located sufficiently deep inside the bubble for the Higgs field to have a constant value. At 
late times, the wall moves with approximately the speed of light. It is expected that the wall 
trajectory and the lower side of the FRW ‘triangle’ will converge asymptotically as the AdS 
boundary is approached, as depicted in fig. 19. The singularity developing below the Cauchy 
horizon appears on a spacelike curve emanating from the point at which the wall reaches the AdS 
boundary. The homogeneity of space viewed by the FRW observer indicates that this singularity 
should correspond to a constant-t surface. In fig. 19 we depict the slicing of the FRW ‘triangle’ 
with such surfaces. The thick solid line represents the possible location of the ‘crunch’, very 
close to the upper side of the FRW ‘triangle’. 

The most important consequence of the above picture is that the ‘crunch’ never reaches the 
bubble wall. This is apparent in fig. 19, as the black solid line, representing the wall trajectory, 
and the blue solid line, representing the ‘crunch’, never cross. They seem to merge on the AdS 
boundary. However, this is an illusion created by the Penrose diagram. The bubble wall lies 



Figure 19: The AdS interior of 
the bubble in conformal coordinates, 
showing the crunch and (in color) a 
patch in FRW coordinates. 
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always slightly outside the FRW ‘triangle’, as its speed never becomes exactly equal to that of 
light. 

A hnal observation relevant for the asymptotic wall expansion concerns the corresponding 
time scales in the various frames. Let us consider a very large bubble with ei = €2 = 1, see 
eqs. (80)-(81), expanding almost at the speed of light, such that R S> £in,2GM and i? S> 1. 
The evolution of the wall in terms of the three time coordinates of the systems (43), (41), (42) 
is given by 

R = =Ro+t, (105) 

1 - C2V 

with Cl ~ 1/(2 £?^ac), C 2 ^ for It is apparent that the wall reaches the AdS 

boundary within a hnite amount of the time coordinate r/, while it requires an inhnite amount 
of time r or t. In particular, an observer in the asymptotically flat spacetime inhnitely far from 
the bubble is reached by the wall only after an inhnite amount of time t. 

A.6 Critical bubbles 

For given and k, corresponding to given interior vacuum energy and surface tension, there 
is a critical bubble radius i?cr- Bubbles that start with negligible wall velocity and R > R^ 
follow trajectories with increasing R, while bubbles that start with R < R^x have diminishing 
R and eventually collapse to a black hole. The critical radius corresponds to the maximum of 
the potential of hg. 15. One can imagine a horizontal line, tangent to the top of the potential. 
The right part of the line is the limiting case of lines starting on the potential and describing 
expanding bubbles, while its left part is the limiting case of similar lines describing collapsing 
bubbles. As the ‘energy’ E is negative, while the potential has a maximum with a positive value 
for Cm = — 1 (see the third plot of hg. 15), it is obvious that there are no critical bubbles with 
negative mass. 

The maximum of the potential and its value at this point are given by eqs. (99)-(100). The 
value of the ‘energy’ can be obtained from eq. (91) with dR/df = 0 and R = i?max- Expressions 
for Rcr and the corresponding mass can then be obtained by combining eqs. (92)-(96). As 
these expressions are not very illuminating we do not present them explicitly. The quantities 
Rcr/^mi GMxxjhvL are functions of the dimensionless combination ntm- In hg. 20 we present the 
functions Rcxinim) ! hn-, GMcxintm) / and [GMcx/i\n]{Rcr/^\n)- 

The critical bubbles have certain characteristics: 

• Their radius is always larger than the Schwarzschild radius. This can be deduced from 
hg. 15, in which it is apparent that the location of the maximum of the potential is always 
outside the horizon. 

• There are two branches of critical bubbles, corresponding to 1/^f^ > or e = 1 (solid 

lines), and or e = — 1 (dashed lines). 

• The radius diverges for k —a l/£in, as the effective energy density in the interior of the 
bubble vanishes in this limit. 
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Figure 20: The radius and mass of critical bubbles as a function of for A > 0 (blue 
curves) and A < 0 (red dashed). 

• The branch with e = 1 reproduces correctly for k —?■ 0 the Newtonian limit of nonrela- 

tivistic bubbles with Ra- = and GMcr = (27. 

• The branch with e = — 1 is not visible to an observer located in region I of the Penrose 
diagram. 

• For given vacuum-energy scale 4n and critical-bubble radius Rcr, the bubbles with e = — 1 
are more massive than the ones with e = 1. (Note that the two types of bubbles also have 
different surface tension k.) 

The most interesting solutions are those that describe bubbles visible to an observer in the 
asymptotically flat region. These are bubbles for which a Newtonian limit exists within their 
parameter range. Their mass-to-radius relation is depicted in the third plot of fig. 20. The 
critical bubbles correspond to the solid line. The parameter range above this line corresponds to 
collapsing bubbles, while the range below to expanding bubbles. Expanding bubbles can have 
negative mass, so their parameter range includes the region below the positive i?-axis. 

A.7 Bubbles in asymptotically de Sitter spacetime 

The evolution of an AdS bubble within an asymptotically dS spacetime can be analysed in 
complete analogy to the previous discussion for an asymptotically flat spacetime. The metric 
of eq. (42) now contains the function font{r) = 1 — ''’^/•^out ~ 2GM/r. There are two horizons, 
corresponding to the zeros of foutir). The Penrose diagram of the Schwarzschild-de Sitter (SdS) 
spacetime is depicted in the right part of fig. 7. It is a combination of the diagrams for the 
Schwarzschild and dS spacetimes [9]. Thick blue lines denote curvature singularities, the dashed 
lines horizons and the dotted lines conformal infinities. The two thin vertical lines at the ends 
of the diagram indicate that the pattern is repeated indefinitely on either side. 

The matching across the domain wall, located at R{t), proceeds as before. We do not analyse 
the many possible cases, as the analysis is a straightforward generalisation of the discussion in 
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the previous subsections. We focus instead on the novel aspects of the SdS case. For a positive 
mass M, the motion of the wall is again determined by eq. (91), with R = pR. However, the 
‘potential’ now has the form 


V{R) 


(t)’ 


^ - 5“^ R^, 
R 


(106) 


with 

= 

where 

A = + j2 K . (108) 

^in ^out 

The form of the ‘potential’ is very similar to that in fig. 15. The horizon corresponds to a value 
Rh such that 

E = -4 - 6^Rjr. (109) 

Rh 

It is again determined by eq. (102), but now has a different shape. In fig. 21 we depict the 
‘potential’ (solid black line) and the horizon (dashed blue line) for 7 = 3, 5 = 1, e = 1. 

The various trajectories correspond to so¬ 
lutions of constant E = — /{G'^M‘^p'^), as 

depicted in fig. 21. The two types of be¬ 
haviour, characterised by e = ± 1 , correspond 
now to the gravitational self-energy of the 
bubble being smaller or larger than the to¬ 
tal difference in energy density 1 /.^?^ -|- 1 / 
between the dS and AdS spacetimes. We do 
not analyse the form of all the possible trajec¬ 
tories, as they are similar to those discussed 
earlier. In fig. 21 we have depicted a few char¬ 
acteristic cases for e = 1 . 

Line A corresponds to a bubble that starts 
below the inner horizon, crosses it, reaches a 
maximal radius and collapses falling again be¬ 
hind the horizon. Figure 21: The ‘potential’ of eq. (106) with 

Line C corresponds to a large bubble that 7 = 3,(5 = l,e = l, for AdS bubbles in 
stars with infinite radius, moves within the Schwarzschild-de Sitter space. 
outer horizon, reaches a minimal radius and 

then re-expands moving again outside the outer horizon. One may consider also the scenario in 
which the bubble is spontaneously created with vanishing wall velocity at a certain radius and 
expands, with the wall moving outside the outer horizon. In this scenario, line C is covered only 
once. 

Line B corresponds to a bubble that starts with a very small radius and expands indefinitely, 
with its wall crossing the inner and outer horizons successively. Its speed asymptotically ap- 



R 
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proaches the speed of light. The form of the wall trajectory on the Penrose diagram is depicted 
in fig. 7. The total space is constructed by patching the white regions of the two plots in fig. 7. 

Line D corresponds to the evolution of a bubble that does not cross any horizons. The reason 
is that ‘energies’ that approach zero correspond to increasing values of the mass parameter M. 
For sufficiently large M, the metric function foutii^) does not vanish at any r, but stays always 
negative. The space has a naked spacelike singularity at r = 0. However, this part of spacetime 
is eliminated and replaced by the interior of the AdS bubble. 

There are many other possibilities for e = — 1 or for negative bubble mass. These can be 
analysed in complete analogy to the trajectories depicted in the second and third plot of fig. 15. 
They correspond to collapsing bubbles or bubbles expanding behind horizons, which are not 
visible to an observer located in the asymptotic de Sitter space. 

The crucial question pertinent to the scenario of Higgs fluctuations during inflation is whether 
the expanding AdS bubbles can completely eliminate the surrounding dS space and thus termi¬ 
nate inflation. It is apparent from fig. 7 that asymptotically the wall trajectory reaches spacelike 
infinity. The wall location separates two spacelike regions: one of them is replaced by the interior 
of the AdS bubble, while the other remains part of an external dS spacetime. Asymptotically 
the total spacetime contains large AdS bubbles within large dS regions. This scenario is in con¬ 
trast to the case of asymptotically flat spacetime, in which the wall asymptotically reaches null 
infinity and the whole space is engulfed by the AdS bubbles. In other words, the inflationary 
growth guarantees that, even when the size of the AdS regions grows with the speed of light, 
the external regions grow even faster, so that they survive at late times. 
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